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ABSTRACT 

Correlations between the intrinsic shapes of galaxies and the large-scale galaxy density field provide an important tool to investigate 
galaxy intrinsic alignments, which constitute the major potential astrophysical systematic in cosmological weak lensing (cosmic 
shear) surveys, but also yield insight into the formation and evolution of galaxies. We measure galaxy position-shape correlations in 
the MegaZ-LRG sample for more than 800,000 luminous red galaxies for comoving transverse separations of 0.3 < r p < 60/r'Mpc, 
making the first such measurement with a photometric redshift sample. In combination with a re-analysis of several spectroscopic 
SDSS samples, we constrain an intrinsic alignment model for early-type galaxies over long baselines in redshift (z < 0.7) and 
luminosity (4 mag) with high statistical precision. 

We develop and test the formalism to incorporate photometric redshift scatter in the modelling of these observations. For r p > 
6/r'Mpc, the fits to galaxy position-shape correlation functions are consistent with the scaling with r p and redshift of a revised, 
nonlinear version of the linear alignment model (Hirata & Seljak 2004) for all samples. An extra redshift dependence oc (1 + j)**" is 
constrained to pother = -0.3 ± 0.8 (lcr). To obtain consistent amplitudes for all data, an additional dependence on galaxy luminosity 
oc If with p = l.Lgi ' s required. The normalisation of the intrinsic alignment power spectrum is found to be (0.077 ± 0.008) p^ 1 for 
galaxies at redshift 0.3 and r band magnitude of -22 (k- and evolution-corrected to z = 0). 

Assuming zero intrinsic alignments for blue galaxies, we assess the bias on cosmological parameters for a tomographic CFHTLS- 
like lensing survey given our new constraints on the intrinsic alignment model parameter space. Both the resulting mean bias and 
its uncertainty are smaller than the lcr statistical errors when using the constraints from all samples combined. The addition of 
MegaZ-LRG data is critical to achieving constraints this strong, reducing the uncertainty in intrinsic alignment bias on cosmological 
parameters by factors of three to seven. 

Key words, cosmology: observations - gravitational lensing: weak - large-scale structure of Universe - cosmological parameters - 
galaxies: evolution 



1. Introduction 

Within the past decade, cosmic shear, the weak gravitational 
lensing effect by the large-scale structure of the Universe, has 
evolved from its first detections (Bacon et al. 2000; Kaiser et al. 
2000; van Waerbeke et al. 2000; Wittman et al. 2000) into a well- 
established and regularly used cosmological probe (see for in- 
stance Benjamin et al. 2007; Fu et al. 2008; Schrabback et al. 
2010 for recent results). Cosmic shear is considered one of the 
most promising techniques to unravel the properties of dark mat- 
ter, dark energy and new gravitational physics (e.g. Albrecht 
et al. 2006; Peacock et al. 2006; Kitching et al. 2008; Thomas 
et al. 2009; Daniel et al. 2010) and acts as a major science driver 
for ongoing, upcoming and future ambitious surveys like Pan- 
STARKS 1 , DES 2 , LSST 3 , JDEM 4 , and Euclid 5 . 



1 Panoramic Survey Telescope & Rapid Response System, 
http : / /pan-starrs . if a . hawaii . edu 

2 Dark Energy Survey, https : //www . darkenergysurvey . org 

3 Large Synoptic Survey Telescope, http : / /www . lsst . org 

4 Joint Dark Energy Mission, http : // j dem . gdf c . nasa . gov 

5 http : //sci . esa. int/science-e/www/area/ 
index. cfm?fareaid=102 



The strength of weak lensing lies in its sensitivity to both the 
geometry of the Universe and the growth of structure, as well as 
its 'clean' relation to the underlying theory which is, except for 
the smallest scales, solely determined by gravitational interac- 
tion. Hence, most efforts in the past years to prepare the analy- 
sis of cosmic shear data for high-precision measurements have 
concentrated on technical issues, in particular the unbiased ex- 
traction of shapes from galaxy images (e.g. Bridle et al. 2010; 
Bernstein 2010; Cypriano et al. 2010; Voigt & Bridle 2010) and 
the estimation of galaxy redshifts from multi-colour photome- 
try (e.g. Abdalla et al. 2007, 2008; Zhang et al. 2010; Bernstein 
& Huterer 2010). Most astrophysical sources of systematic er- 
rors such as the limited accuracy of N-body simulations (e.g. 
Hilbert et al. 2009) and baryonic effects on structure growth (e.g. 
White 2004) are currently subdominant and can, like the techni- 
cal systematics, largely be reduced by more powerful simula- 
tions, specifically designed instrumentation, and optimised anal- 
ysis techniques. 

However, this statement does not apply to the intrinsic align- 
ment of galaxies, which is a major astrophysical systematic af- 
fecting all two-point and higher-order weak lensing statistics. 
The estimation of gravitational shear correlations from the ob- 
served ellipticities of galaxy images is simple if one assumes 
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that the intrinsic ellipticity of one galaxy is not correlated with 
either the intrinsic ellipticity of or the gravitational shear act- 
ing on another galaxy. Both of these assumptions are incorrect 
in the presence of intrinsic alignments which are e.g. caused by 
tidal torquing and stretching of galaxies by the large-scale grav- 
itational field. These effects depend intricately on galaxy for- 
mation and evolution including baryonic physics, thus hamper- 
ing modelling via analytical calculations or N-body simulations. 
The bias on cosmological parameters when ignoring intrinsic 
alignments can be severe; e.g. for the dark energy equation of 
state parameter, Bridle & King (2007) find a deviation of 50 % 
from a fiducial value of -1 for a Euclid-like survey if intrinsic 
alignment contamination is not accounted for in the analysis. 

Correlations between the intrinsic ellipticities of two galax- 
ies (hereafter II) can occur if the galaxies are subject to the tidal 
forces of the same local or large-scale matter structures. Thus, 
a pair of galaxies with II correlations must be physically close, 
i.e. at small separation both in redshift and on the sky, which 
allows for a relatively straightforward removal of the II signal 
from cosmic shear data if photometric redshifts are sufficiently 
accurate (King & Schneider 2002, 2003; Heymans & Heavens 
2003; Takada & White 2004). The mutual alignment of halo 
shapes and spins has been studied extensively in N-body simula- 
tions (Splinter et al. 1997; Onuora & Thomas 2000; Faltenbacher 
et al. 2002; Hopkins et al. 2005; Faltenbacher et al. 2007, 2008; 
Lee et al. 2008), to a limited degree also including the effect of 
baryonic physics (van den Bosch et al. 2002; Bett et al. 2010; 
Hahn et al. 2010). Models for II correlations have been devel- 
oped analytically or via fits to simulations (Croft & Metzler 
2000; Heavens et al. 2000; Catelan et al. 2001; Crittenden et al. 
2001; Jing 2002; Mackey et al. 2002; Hirata & Seljak 2004; 
Schneider & Bridle 2010), with widely varying results. The II 
signal can either be observed using galaxy ellipticity correlations 
at very low redshift where cosmic shear is negligible (Brown 
et al. 2002 in SuperCOSMOS) or by selecting galaxy pairs with 
small physical separation (Mandelbaum et al. 2006; Hirata et al. 
2007; Okumura et al. 2009; Brainerd et al. 2009 in various SDSS 
samples). 

Hirata & Seljak (2004) identified gravitational shear- 
intrinsic ellipticity cross-correlations (hereafter GI) as a further, 
potentially more serious contaminant of cosmic shear surveys. 
GI correlations are generated by matter that aligns a nearby fore- 
ground galaxy and at the same time contributes to the gravita- 
tional lensing signal of a background galaxy. Thus, GI correla- 
tions are not restricted to close (along the line-of-sight) pairs. 
This type of intrinsic alignment has a redshift dependence that is 
very similar to that of lensing, so that GI correlations could be 
particularly prominent in the deep data of future cosmic shear 
surveys. The underlying correlations between halo shapes or 
orientations and the surrounding matter structure have been de- 
tected on a large range of mass scales in simulations (Bailin & 
Steinmetz 2005; Altay et al. 2006; Heymans et al. 2006; Kuhlen 
et al. 2007). However, as for II correlations, modelling attempts 
(Hui & Zhang 2002; Hirata & Seljak 2004; Heymans et al. 2006; 
Schneider & Bridle 2010) do not currently yield robust predic- 
tions. 

In the absence of a compelling model for GI correlations, 
methods with a weak dependence on the intrinsic alignment 
model are required to eliminate biases on cosmology in cosmic 
shear analyses. Joachimi & Schneider (2008, 2009) made use of 
the characteristic redshift dependence of the GI signal to null it 
in a fully model-independent way, while King (2005), Bridle & 
King (2007), Bernstein (2009), and Joachimi & Bridle (2010) in- 
troduced very general parametrisations of the intrinsic alignment 



contributions containing parameters that are then marginalised 
over (see Kirk et al. 2010 for an application to data). Both ap- 
proaches feature a considerable loss of information, substan- 
tially weakening constraints on cosmological model parameters. 
Observational constraints on the GI contribution are thus crucial, 
as they can tighten priors on the intrinsic alignment parameters 
that need to be marginalised over, and shed light on the underly- 
ing physical processes by constraining intrinsic alignment mod- 
els. 

To study the GI signal, one must measure the correlations 
between the matter distribution and the intrinsic shear, i.e. the 
correlated part of the intrinsic galaxy shapes (see e.g. Hirata & 
Seljak 2004). Assuming linear biasing, and constructing corre- 
lation functions with negligible gravitational lensing contribu- 
tions, this measurement involves cross-correlating galaxy num- 
ber densities and ellipticities (Mandelbaum et al. 2006; Hirata 
et al. 2007; Okumura & Jing 2009; Mandelbaum et al. 2010). 
These galaxy number density-intrinsic shear correlations were 
detected by Mandelbaum et al. (2006) in SDSS spectroscopic 
data at low redshift (z ~ 0.1). Hirata et al. (2007) extended the 
analysis to SDSS Luminous Red Galaxies (LRGs) at slightly 
higher redshift (0.15 < z < 0.35), finding evidence for a depen- 
dence on galaxy type and for an increase of the intrinsic align- 
ment signal of early-type galaxies with luminosity. They also 
considered a small sample of galaxies at intermediate redshifts 
(z ~ 0.5) from the 2dF-SDSS LRG and Quasar Survey (2SLAQ, 
Cannon et al. 2006), but only marginally detected a signal in a 
bright subsample due to poor statistics. Recently, Mandelbaum 
et al. (2010) increased the redshift baseline for position-shape 
correlation measurements of blue galaxies out to z ~ 0.7 by in- 
corporating a large number of spectroscopic redshifts from the 
WiggleZ Survey (Drinkwater et al. 2010), reporting a null detec- 
tion for all subsamples. 

Analogous to Mandelbaum et al. (2010), this work presents 
observational constraints for early-type galaxies (out to z < 0.7), 
using SDSS shape measurements together with photometric red- 
shift information from the MegaZ-LRG catalogue (Collister 
et al. 2007; Abdalla et al. 2008). In combination with previously 
analysed red galaxy samples, we cover a wide range of redshifts 
and galaxy luminosities with high statistical precision, which al- 
lows us to narrow down the redshift and luminosity evolution of 
galaxy number density-intrinsic shear correlations. The longer 
baselines in redshift and luminosity enable a meaningful extrap- 
olation to typical parameter values found for galaxies in cosmic 
shear surveys, so that we can estimate the contamination due to 
intrinsic alignments for present-day surveys. For the first time, 
we include a galaxy sample that only has photometric redshift 
information into the analysis, and we develop the corresponding 
formalism. In particular, we account for the spread of the num- 
ber density-intrinsic shear correlations along the line of sight due 
to photometric redshift scatter, and determine the importance of 
other signals such as galaxy-galaxy lensing in the observables. 

This paper is structured as follows: In Sect. 2 we provide an 
overview of the galaxy samples used in the analysis. Section 3 
contains the methodology for the correlation function measure- 
ment. We develop the modelling of photometric redshift number 
density-intrinsic shear correlations in Sect. 4, deferring the tech- 
nical aspects of the derivation, as well as a revision of the red- 
shift dependence of the linear alignment model, to two appen- 
dices. In Sect. 5 the results of our analysis are presented, includ- 
ing constraints on an intrinsic alignment model and a discussion 
of systematic tests. These findings are then applied in Sect. 6 to a 
prediction of the intrinsic alignment contamination of a generic 
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present-day tomographic cosmic shear survey, before we sum- 
marise and conclude in Sect. 7. 

Where appropriate, we follow the notation of Mandelbaum 
et al. (2006), Hirata et al. (2007), and Mandelbaum et al. (2010). 
Since in contrast to these foregoing works we have to con- 
sider various contributing signals, we denote the observable by 
'galaxy position-shape' correlations and otherwise follow the in- 
dexing scheme of Joachimi & Bridle (2010). Throughout we will 
assume as our cosmological model a spatially fiat ACDM uni- 
verse with parameters Q m = 0.25, Qb = 0.05, erg = 0.8, h = 0.7, 
and n s = 1 .0. While h = 0.7 is used for e.g. power spectrum cal- 
culations, all distances etc. will be given in units that are inde- 
pendent of h, i.e. in /z^'Mpc. We use the AB magnitude system 
and specify luminosities of the galaxy samples under consider- 
ation in the SDSS r filter. Absolute magnitudes are consistently 
given in terms of h = 1 and typically k + e-corrected to z — 
unless stated otherwise. 

2. Data 

All data used in this paper come from the Sloan Digital Sky 
Survey (SDSS). The SDSS (York et al. 2000) imaged roughly 
n steradians of the sky, and followed up approximately one mil- 
lion of the detected objects spectroscopically (Eisenstein et al. 
2001; Richards et al. 2002; Strauss et al. 2002). The imaging 
was carried out by drift-scanning the sky in photometric condi- 
tions (Hogg et al. 2001; Ivezic et al. 2004), in five bands (ugriz) 
(Fukugita et al. 1996; Smith et al. 2002) using a specially- 
designed wide-field camera (Gunn et al. 1998). These imaging 
data were used to create the galaxy shape measurements that we 
employ in this paper. All of the data were processed by com- 
pletely automated pipelines that detect and measure photomet- 
ric properties of objects, and astrometrically calibrate the data 
(Lupton et al. 2001; Tucker et al. 2006). The SDSS has had 
seven major data releases, and is now complete (Stoughton et al. 
2002; Abazajian et al. 2003, 2004, 2005; Finkbeiner et al. 2004; 
Adelman-McCarthy et al. 2006, 2007, 2008; Abazajian et al. 
2009). 

2.1. MegaZ-LRG 

The MegaZ-LRG sample (Collister et al. 2007) is based on 
SDSS five-band (ugriz) imaging data. It contains more than 
a million luminous red galaxies at intermediate redshifts be- 
tween 0.4 and 0.7, i.e. beyond the redshifts of the LRGs al- 
ready targeted with the SDSS spectrograph (z < 0.45, Eisenstein 
et al. 2001). While the original catalogue was selected from the 
4th SDSS data release, we use the updated version presented 
in Abdalla et al. (2008) based on data release 6 (Adelman- 
McCarthy et al. 2008). 

The photometry in five bands is used to determine photo- 
metric redshifts for the MegaZ-LRG sample. For a subset of the 
galaxies, spectroscopic redshift information is required for cal- 
ibration and cross-checking, which is provided by the 2SLAQ 
survey (Cannon et al. 2006). Consequently, the selection criteria 
of MegaZ-LRG have been designed to match those of 2SLAQ, 
using a series of magnitude and colour cuts (for details see 
Cannon et al. 2006; Collister et al. 2007). These criteria have 
an efficiency of 95 % in detecting LRGs in the redshift range 
0.4 < z < 0.7, the failures being almost entirely due to M-type 
stars. In Sect. 3 we will describe how we account for this con- 
tamination in our analysis. 

The 2SLAQ selection criteria fluctuated a little at the be- 
ginning of the survey. Specifically, the faint limit of the i 
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Fig. 1. Top panel: Redshift distributions of the galaxy samples analysed. 
Shown are the histograms for the SDSS LRG samples in black (faint, 
M r > -22.3), blue (medium, -22.6 < M r < -22.3), and purple (bright, 
M r < -22.6), for the full MegaZ-LRG sample in red, and for the SDSS 
Main samples in green (L4) and orange (L3). Note that both SDSS LRG 
and MegaZ-LRG samples are split into two redshift bins each, the SDSS 
LRG samples at z = 0.27 and the MegaZ-LRG sample at z = 0.529. 
Bottom panel: Distribution of rest-frame absolute magnitudes M r . The 
colour coding of the histograms is the same as in the top panel. In the 
case of the SDSS LRG and MegaZ-LRG samples solid lines correspond 
to the low redshift bin, dotted lines to the high redshift bin, respectively. 
Note that the MegaZ-LRG histograms rely on photometric redshift esti- 
mates, and that they have been downscaled by a factor of 20 to facilitate 
the comparison with the other samples. 



band magnitude ;'d e v (the total magnitude estimated using a de 
Vaucouleurs profile), and the minimum value of d ± — (r - i) - 
(g ~ r )/8 (a colour variable used to select LRGs), were var- 
ied slightly. For the majority of the 2SLAQ survey, the crite- 
ria /dev ^ 19.8 and d ± > 0.55 were used. For further details 
on this see Cannon et al. (2006). However, for the full MegaZ- 
LRG sample described in Collister et al. (2007), the flux limit 
is i'dev ^ 20, which means that roughly 1/3 of the sample is 
fainter than the 2SLAQ flux limit. Details about the photometric 
redshift estimation are provided in Sect. 2.4. 

As summarised in Table 1, we use about 860,000 galaxies 
with a mean redshift of 0.54 in the full MegaZ-LRG sample to 
compute galaxy number densities. The total number of galaxies 
is less than that of the full MegaZ-LRG catalogue by Collister 
et al. (2007) because a fraction of the area comes from imag- 
ing data that was not yet processed by the shape measurement 
pipeline used to estimate galaxy shapes. As will be discussed 
in Sect. 2.5, accurate shape measurements could be obtained for 
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sample 
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212944 


U.J7 


1 OS 


SDSS LRG, faint, z < 0.27 


4314 


16701 


7030 


0.21 


1.06 


SDSS LRG, faint, z > 0.27 




19397 


9038 


0.32 


1.07 


SDSS LRG, medium, z < 0.27 




16701 


6139 


0.22 


1.50 


SDSS LRG, medium, z > 0.27 




19397 


6700 


0.31 


1.50 


SDSS LRG, bright, z < 0.27 




16701 


3532 


0.22 


2.13 


SDSS LRG, bright, z > 0.27 




19397 


3659 


0.31 


2.12 


SDSS Main L4 red 


3423 


288863 


26872 


0.10 


0.33 


SDSS Main L3 red 




288863 


10867 


0.07 


0.14 



Notes. The columns give the sample name, the survey area covered, the number of galaxies Af ga i used to compute both the galaxy number densities 
and the shapes, the mean redshift (z), and the mean luminosity (L) /L in terms of a fiducial luminosity Lq that corresponds to a k + e-corrected (to 
z = 0) absolute r band magnitude M = -22. 



about 50 % of the MegaZ-LRG galaxies, which are then used to 
trace the intrinsic shear field. 

We also compute r band luminosities, taking into account 
galactic dust extinction (using reddening maps from Schlegel 
et al. 1998 and extinction-to-reddening ratios from Stoughton 
et al. 2002), and £+e-correcting the model magnitudes to z = 
using the same templates as in Wake et al. (2006). Luminosities 
are given in terms of a fiducial luminosity Lq, corresponding to a 
fc+e-corrected absolute r band magnitude of Mo = -22 mag. Due 
to this procedure, the corrected luminosity acts as a tracer for 
the total stellar mass, and consequently also for the total mass, 
of the galaxy. For reference, Blanton et al. (2003b) obtain a rest- 
frame magnitude at z — 0. 1 of M* = -20.44 for a purely flux- 
limited sample of SDSS spectroscopic galaxies at z < 0.2. Using 
Wake et al. (2006) to k + e-correct this value to z = 0, we find 
M* r = -20.32 or L* = 0.21L () . 

Due to a magnitude-dependent redshift success rate and a 0.2 
magnitude difference in the i^v cuts (see above), the MegaZ- 
LRG sample is somewhat fainter in absolute magnitude than 
the 2SLAQ sample. Thus, we define the z = absolute mag- 
nitudes using the MegaZ-LRG photometric redshift values, but 
we also multiply them by a correction factor that is derived from 
2SLAQ. In 2SLAQ, we find that if we use the photometric red- 
shift to define the luminosity rather than the spectroscopic red- 
shift, the mean sample luminosity appears to be 5% too low. 
Possible sources of this bias are discussed in Sect. 2.4. Thus, we 
correct the mean (L) /Lq for the MegaZ-LRG sample, as defined 
using photometric redshifts, upwards by 5% to account for the 
impact of photometric redshift errors, the corrected values being 
given in Table 1 . 

In addition, with a cut at photometric redshift z = 0.529, we 
split the sample into two redshift bins, each containing roughly 
the same number of galaxies. We show the MegaZ-LRG photo- 
metric redshift distribution, as well as the luminosity distribu- 
tions of both MegaZ-LRG subsamples in Fig. 1, where the latter 
are also based on photometric redshift estimates. As is evident 
from the figure, the redshift cut for the MegaZ-LRG sample also 
segregates the galaxies in luminosity since the sample is magni- 
tude limited. 

2.2. Spectroscopic LRGs 

We also use the SDSS DR4 spectroscopic LRG sample 
(Eisenstein et al. 2001), for which measurements of galaxy 



position-shape correlations were presented in Hirata et al. 
(2007). No new measurements of this sample are made for the 
current work; instead, we combine the previous measurements 
with the new MegaZ-LRG data to provide constraints over a 
longer redshift baseline, see Fig. 1 and Table 1 . 

The SDSS spectroscopic LRG sample has a flux limit of r < 
19.2 and colour cuts to isolate LRGs. We include these galaxies 
in the redshift range 0.16 < z < 0.35, for which the sample 
is approximately volume-limited, and includes 36 278 galaxies 
total. 

In order to study variation within this sample, we use cuts on 
several parameters. First, we construct luminosities using the r 
band model magnitudes in the same way as for the MegaZ-LRG 
sample (Sect. 2.1), and define three luminosity subsamples with 
M r < -22.6 ('bright'), -22.6 < M r < -22.3 ('medium'), and 
M r > -22.3 ('faint'), see Table 1. The absolute magnitude cuts 
are defined in terms of h — Ho /(100 km s _1 Mpc~') such that one 
can implement the cuts without specifying the value of Hq. The 
magnitudes have been corrected for galactic extinction, and are 
k + e-corrected to z = using the same templates as in Wake 
et al. (2006). Each luminosity subsample is in addition split at 
z = 0.27 into a high- and a low-redshift bin. 

2.3. Main sample 

Furthermore we consider galaxies at a more typical luminos- 
ity, the same subsamples of the DR4 SDSS Main spectroscopic 
sample as in Mandelbaum et al. (2006), divided by luminos- 
ity and other properties. For this work, we use the red galaxies 
in L3 (from roughly 1 .5 to 0.5 magnitudes fainter than typical 
L* galaxies) and L4 (from 0.5 magnitudes fainter to 0.5 magni- 
tudes brighter than L,) 6 . The sample properties were described 
in full in that paper; for this work, we mention only that the 
luminosities described to select the sample are Petrosian magni- 
tudes, extinction-corrected using reddening maps from Schlegel 
et al. (1998) with the extinction-to-reddening ratios given in 
Stoughton et al. (2002), and fc-corrected to z = 0.1 using kcor- 
rect v3_2 software as described by Blanton et al. (2003a). For 
consistency with previous work, L3 and L4 were initially se- 
lected with respect to this type of absolute magnitude. However, 
in order to compare this sample with the others, we must also 



6 We refrain from including the L5 and L6 samples as well because 
of their significant overlap with the spectroscopic LRG sample. 
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compute absolute magnitudes in the same way as for those 
(model instead of Petrosian magnitudes, k + e-corrected to z = 
instead of 0.1 using specifically red galaxy templates). To use 
the k + e-corrections as in the previous sections, we must first 
ensure that we are selecting a properly red galaxy sample for 
which they are appropriate. 

Hirata et al. (2007) used an empirically-determined redshift- 
dependent colour separator of u - r = 2.1 + 4.2z, which uses 
observer-frame rather than rest-frame colours; within these lu- 
minosity bins, the fractions of red galaxies were 0.40 (L3), 0.52 
(L4), 0.64 (L5), and 0.80 (L6). For the following analysis, we 
want to define a red-sequence sample that is comparable to the 
higher-redshift samples defined previously (which are typically 
more strictly defined, though see Wake et al. 2006 and Sect. 5.4 
for a discussion of the differences between MegaZ-LRG and 
SDSS LRGs). Thus, for this work we use a different colour sep- 
arator and recompute the GI correlations for the new 'red' sam- 
ples. 

To generate the new colour separator, we first fc-correct the 
extinction-corrected model magnitudes using kcorrect to z — 0. 
The reason for beginning in this fashion, rather than using the 
Wake et al. (2006) templates as for the SDSS LRG and MegaZ- 
LRG samples, is that those k + e-corrections are not actually 
applicable to the majority of the Main sample since it is flux- 
limited and therefore has a significant fraction of blue galax- 
ies. Thus, we begin with kcorrect which imposes a template- 
dependent ^-correction. Then, we place a cut in the distribution 
of rest-frame g — r, to strictly isolate the reddest galaxies. The 
cut value was motivated by Padmanabhan et al. (2004), which 
shows typical values for pressure-supported elliptical galaxies in 
the Main sample; however the actual cut value is shifted because 
the k-correction in that paper is to z = 0. 1 rather than 0. The new 
red fractions in L3 and L4 are 22 % and 26 %, respectively. For 
the galaxies satisfying this cut, we then get absolute magnitudes 
using the Wake et al. (2006) k + e-corrections, which we find 
are typically 0.2 magnitudes brighter than the original Petrosian 
z — 0.1 magnitudes used to define the L3 and L4 samples. 

The detailed issue of consistency between the different sam- 
ples in this paper (MegaZ-LRG, spectroscopic LRGs, and Main 
L3 and L4 red samples) is irrelevant to any attempts to simply 
fit each sample to an intrinsic alignment amplitude. However, 
before attempting to combine the measurements in the different 
samples, we must address this issue. Thus, a detailed comparison 
will be given in Sect. 5.4. 

2.4. Photometric redshifts 

The MegaZ-LRG catalogue is selected from the SDSS imaging 
database using a series of colour and magnitude cuts (Collister 
et al. 2007) that match the default selection criteria of the 2SLAQ 
survey (Cannon et al. 2006) except for going 0.2 magnitudes 
fainter. The spectroscopic redshifts available from 2SLAQ were 
used to train and test the photometric redshift code, which was 
applied to the entire set of LRGs selected from the SDSS imag- 
ing database. 

Around 13,000 spectroscopic objects were obtained in the 
2SLAQ survey, of which about 8000 were used to train a neu- 
ral network (Collister & Lahav 2004) to produce photometric 
redshifts, leaving approximately 5100 galaxies to verify the esti- 
mates (Collister et al. 2007; Abdalla et al. 2008). The quality of 
the photometric redshifts is very good because of the large 4000 
Angstrom break present in luminous red galaxies. Accurate pho- 
tometric redshifts are paramount for our analysis in order to ob- 
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Fig. 2. Binned histogram of spectroscopic redshifts from 2SLAQ and 
photometric redshift estimates from the MegaZ-LRG catalogue. Note 
that the shading of the bins is logarithmic. The solid line indicates a 
one-to-one relation between spectroscopic and photometric redshifts, 
coinciding with the mean trend to high accuracy. The dotted lines cor- 
respond to the ± l<r scatter. 



tain good signal-to-noise, and to ensure that the matter-intrinsic 
correlations dominate the observed signals, see Sect. 4.2. 

We plot in Fig. 2 the quality of the photometric redshifts that 
can be tested with the verification sample of ~ 5000 2SLAQ 
spectra. We find that the distribution of the differences between 
photometric redshift estimate and spectroscopic redshift has a 
mean of zero. As a consequence, the mean trend of the photo- 
metric redshift distribution, given a spectroscopic redshift, is in- 
distinguishable from a one-to-one relation. The number of galax- 
ies with a difference between photometric redshift I and spec- 
troscopic redshift z largely exceeding the typical scatter (5cr or 
more) is always less than 3 % for any photometric redshift bin in 
the range z < 0.65. The distribution of differences between pho- 
tometric and spectroscopic redshift is well fit by a Gaussian with 
width 0.024(1 +z) (corresponding to the scatter shown in the fig- 
ure), in good agreement with the results by Collister et al. (2007) 
who find a very similar scatter in the range 0.45 < z < 0.50 
in which most galaxies of our sample reside. Their scatter in- 
creases by up to 50 % for higher photometric redshifts in the 
range 0.60 < z < 0.65. 

While the distribution of spectroscopic redshifts given a pho- 
tometric redshift, i.e. p(z\z) is unbiased on average, the distri- 
butions p(z\z) have significant systematic offsets, as is evident 
from Fig. 2. Photometric redshift estimates for galaxies at low 
redshift are biased high, those for galaxies at high redshift are 
biased low, leading to a distribution of photometric redshifts that 
is more compact than the corresponding spectroscopic distribu- 
tion. Since faint galaxies are preferentially found at low redshift, 
and bright galaxies at high redshift, the luminosity distribution 
of MegaZ-LRG galaxies based on photometric redshifts as dis- 
played in Fig. 1 also appears more compact than the true distri- 
bution. We are not directly affected by this change in the shape 
of the distribution because we base our analysis on mean sam- 
ple luminosities. However, as was discussed in Sect. 2.1, the re- 
distribution of galaxies due to photometric redshifts also leads 
to a small change in the mean luminosity of the MegaZ-LRG 
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Fig. 3. Fraction of galaxies in the MegaZ-LRG sample with high-quality 
shape measurements, as a function of photometric redshift (top left 
panel), rest-frame absolute r band magnitude M r (top right panel), i 
band de Vaucouleurs magnitude ! dcV which was used as a selection cri- 
terion for the MegaZ-LRG catalogue (bottom left panel), and apparent 
observer-frame r band magnitude used to impose a magnitude limit on 
the shape catalogue (bottom right panel). The red histograms show the 
match fraction for the full MegaZ-LRG shape sample, and the blue his- 
tograms for shape sample with the additional colour cut that will be 
discussed in Sect. 5.4. For reference we have added to each panel the 
histogram of the full MegaZ-LRG sample with arbitrary normalisation 
as black dotted lines. Note that the fraction of galaxies with shape in- 
formation does not depend strongly on redshift and r band magnitude. 



samples which we correct for via the photometric versus spec- 
troscopic redshift relation from 2SLAQ. 

Note that we employ the photometric versus spectroscopic 
redshift histogram of Fig. 2 directly in the correlation function 
models as an approximation to p(z\z), despite the 0.2 magnitudes 
fainter catalogue used in this work. We will show in Sect. 5.1 that 
this is not a significant cause of systematic error in our analysis, 
as expected because the photometric redshift properties are only 
slightly extrapolated to /d e v = 20. We note furthermore that, due 
to the small field of the 2SLAQ survey, cosmic variance limits 
the accuracy of using this photometric versus spectroscopic red- 
shift relation for MegaZ-LRG as well. 

2.5. Galaxy shape measurements 

To measure the intrinsic shears, we use PSF-corrected shape 
measurements from SDSS, specifically the galaxy ellipticity 
measurements by Mandelbaum et al. (2005), who obtained 
shapes for more than 30 million galaxies in the SDSS imaging 
data down to extinction-corrected model magnitude r = 2 1 .8 us- 
ing the Reglens pipeline. We refer the interested reader to Hirata 
& Seljak (2003) for an outline of the PSF correction technique 
(re-Gaussianisation) and to Mandelbaum et al. (2005) for all de- 
tails of the shape measurement. The two main criteria for the 
shape measurement to be considered high quality are that galax- 
ies must (a) have extinction-corrected r band model magnitude 
r < 21.8, and (b) be well-resolved compared to the PSF size in 
both r and i bands. 



The fraction of MegaZ-LRG galaxies with high-quality 
shape measurements is 50 %, nearly independent of photomet- 
ric redshift and r band magnitude, see Fig. 3. We also plot this 
match fraction as a function of the i band de Vaucouleurs mag- 
nitude, where i^y < 20 was used as the magnitude cut for the 
MegaZ-LRG catalogue, and do not find a significant evolution 
either. It is interesting to note that the observer-frame r - /d e v 
strongly increases with redshift, exceeding unity around z = 0.5, 
i.e. approximately at the peak of the MegaZ-LRG redshift dis- 
tribution. Therefore all galaxies close to the r band magnitude 
limit have to be located at high redshift, given the limit on i'd e v 
for MegaZ-LRG, so that any evolution of the match fraction with 
redshift should imply an evolution with r. Since both are roughly 
constant in spite of the size cuts in galaxy shape measurements, 
the decrease in apparent galaxy size due to the larger distance at 
higher redshift has to be balanced by an increase of the physical 
dimensions of these galaxies. This can be explained by the pos- 
itive correlation between galaxy size and absolute magnitude, 
where for a flux-limited survey the galaxies at higher redshift 
are intrinsically brighter on average. 

The match fraction as a function of absolute rest-frame r 
band magnitude M r displays a moderate decrease towards intrin- 
sically fainter galaxies. Since in this case the effects of different 
galaxy distances have been removed, the galaxy size-luminosity 
relation causes the more luminous galaxies to have generally a 
higher share of good shape measurements. Note that in Fig. 3 we 
also show match fractions for the MegaZ-LRG sample with an 
additional colour cut that removes the bluest galaxies from the 
sample, see Sect. 5.4. This cut predominantly removes galaxies 
which are faint in the r and i'd e v bands, without a significant de- 
pendence on redshift. 

Moreover we find that the fraction of galaxies with high- 
quality shape measurements exhibits a considerable dependence 
on observing conditions due to the resolution cut. We have gen- 
erated mock catalogues that account for all of these variations in 
the ability to measure galaxy shapes with magnitude and observ- 
ing conditions. 



3. Measurement of correlation functions 

The software for the computation of galaxy position-shape cor- 
relation functions is one of the codes used in Mandelbaum et al. 
(2006), Hirata et al. (2007), and Mandelbaum et al. (2010). 
Here we summarise the methodology, and refer the reader to 
Mandelbaum et al. (2010) for details. This software finds pairs of 
galaxies (with one belonging to the shape-selected sample used 
to trace the intrinsic shear field, and the other belonging to the 
full sample used to trace the density field) using the SDSSpix 
package 7 . We measure the correlations as a function of comov- 
ing transverse separation r p and comoving line-of-sight sepa- 
ration n of the galaxy pairs, over the complete range of red- 
shifts covered by a sample. The correlation functions must be 
computed using a large range of n because of photometric red- 
shift error (see Sect. 4); we divide this range into bins of size 
An = 10/r'Mpc. 

We adopt a variant of the estimator presented in 
Mandelbaum et al. (2006), which is given by 



S+D-S+R 



(1) 



http : //lahmu . phyast . pi tt . edu/" scranton/SDSSPix/ 
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where S +D stands for the correlation between all galaxies in the 
full MegaZ-LRG catalogue, tracing the density field, and those 
from the subset with shape information, 



S+D 



i*j\r p M 



g+010 
2<R 



(2) 



Here, e + (j\i) denotes the radial component of the ellipticity of 
galaxy j measured with respect to the direction towards galaxy 
i out of the number density sample. A similar equation holds 
for S +R, but in this case the galaxies of the number density 
sample are taken from a random catalogue compliant with the 
selection criteria of the MegaZ-LRG data. The subtraction of 
S +R is meant to remove any spurious shear component, as de- 
scribed in previous works. The shear responsitivity % represents 
the response of our particular ellipticity definition to a shear, and 
for an ensemble with rms ellipticity per component of e rms , is 
roughly 1 - e^ ms « 0.85. For more details on the random cat- 
alogue generation and treatment and the shear responsitivity % 
see Mandelbaum et al. (2006). 

The normalisation of (1) is given by the number of pairs D S R 
with one real galaxy from the shape subsample D s and one ran- 
dom galaxy from the full random catalogue R. All these pair 
counts are done for galaxies with transverse comoving separa- 
tion r p and comoving line-of-sight separation n, for all redshifts 
z in the sample. 

To deduce matter-intrinsic shear correlations from £ g+ , we 
also measure galaxy clustering for MegaZ-LRG via 



DsD 
D S R 



1 



(3) 



where DsD denotes the number of galaxy pairs between the full 
MegaZ-LRG catalogue (used to trace the density field) and the 
intrinsic shear sample, and D S R is the number of pairs with one 
real galaxy in the sample used to trace the intrinsic shear and 
one in the (full) random catalogue representing the density field 
sample. Again, all these galaxies are selected from bins in r p 
and n over all redshifts in the sample. By cross-correlating the 
MegaZ-LRG full and shape samples in (3), we intend to miti- 
gate the effect of a residual star contamination /contain = 0.05 in 
the sample (see Sect. 2.1) which should enter (3) only linearly, 
as will be detailed in Sect. 5.2. Determining galaxy clustering 
from auto-correlations of the full galaxy sample would result in 
a contamination which is quadratic in / con tam to leading order. 

The projected correlation function is then computed by sum- 
mation of the correlation functions in (1) and (3) along the line 
of sight, multiplied by All. This calculation is done in Nbi a = 10 
logarithmically spaced bins in transverse separation in the range 
0.3 < r p < 60/z~'Mpc. We re-bin the existing correlation func- 
tions in r p for the SDSS LRG and Main samples accordingly. 
For the spectroscopic samples, bins in line-of-sight separation 
have a width of All = 4/r'Mpc. The cut-off in this stack- 
ing process is Il max = 60/z~ 1 Mpc for the spectroscopic data 
sets, capturing virtually all of the signal (Mandelbaum et al. 
2006; Hirata et al. 2007; Padmanabhan et al. 2007; Mandelbaum 
et al. 2010). MegaZ-LRG correlation functions are computed 
for n max = 90/r'Mpc and n max = 180/i _1 Mpc, where we 
will investigate the effect of these truncations in more detail in 
Sect. 5.1. Note that these values of Il max were chosen to roughly 
coincide with the lcr and 2<x photometric redshift scatter. 

Covariance matrices for MegaZ-LRG are determined using 
a jackknife with 256 regions, in order to account properly for 
shape noise, shape measurement errors, and cosmic variance. 
For the total survey area probed by MegaZ-LRG the maximum 



comoving transverse separations contained in each jackknife re- 
gion are well above 100/r'Mpc at z = 0.5, and thus consider- 
ably larger than the maximum r p used for the analysis, so that the 
jackknife samples are independent. 50 jackknife regions were 
used to obtain covariances for the SDSS LRG and Main samples 
(Hirata et al. 2007), where the smaller number of regions is a 
consequence of the lower mean redshift and the smaller survey 
area covered. We address the issue of noise in the covariance 
matrices below in Sect. 4.4. 

To compute correlation functions as a function of the comov- 
ing separations r p and n, one needs to assume a fiducial cos- 
mology to transform the observable galaxy redshifts and angular 
separations. The cosmological model for this conversion differs 
from our fiducial cosmology in the value Q. m = 0.3 for all sam- 
ples under consideration, in order to maintain consistency with 
the signals computed for the SDSS LRG sample from Hirata 
et al. (2007). Since we use our default value of Q. m = 0.25 for all 
model calculations, this discrepancy could bias our results. We 
evaluate the effect of the difference in O m for the MegaZ-LRG 
high-redshift bin which is the sample at the highest redshift and 
should thus be affected most. The change in n can be safely ne- 
glected while r p changes by 2 % at z = 0.59. Converting the 
transverse separation of the observed correlation functions from 
Q. m = 0.3 to Q m = 0.25, we find changes in the fit results for the 
galaxy bias and the intrinsic alignment model amplitude below 
1 % each, and thus conclude that the discrepancy in £2 m can be 
neglected in our analysis. 

4. Modelling 

While the methodology for the analysis of spectroscopic galaxy 
samples is already well established (Mandelbaum et al. 2006; 
Hirata et al. 2007), we consider for the first time a sample 
with only photometric redshift information, obtained from the 
MegaZ-LRG catalogue. In this section and in Appendix A we 
derive the models that are later compared to the observational 
data, incorporating photometric redshift uncertainty. The pho- 
tometric redshift scatter has two major effects on our measure- 
ments: First, the truncation of the observed correlation function 
at large line-of-sight separations n has to be taken into account 
explicitly in the model. Second, we must assess the importance 
of contributions other than galaxy number density-intrinsic shear 
correlations to the observations. For reasons of optimum signal- 
to-noise and a simplified physical interpretation, the observa- 
tions are expressed as line-of-sight projected correlation func- 
tions as a function of comoving transverse separation between 
galaxy pairs, and we transform the model accordingly. 

4.1. The number density-intrinsic shear correlation function 

As a first step, we compute the correlation function be- 
tween galaxy number density and intrinsic shear (hereafter 
gl), ^ ot (r p ,T\,z m ). As before, f p denotes the comoving trans- 
verse separation, and tl the comoving line-of-sight separation 
of galaxy pairs which are located at a mean redshift z m . Note 
that we assign a bar to all quantities derived from photomet- 
ric redshift estimates. In the following we will refer to correla- 
tions of the form £(r p , II, z) as the three-dimensional correlation 
function. We follow the notation of Joachimi & Bridle (2010) 
in denoting the different signals that contribute to the observa- 
tions 8 . As a reminder, for each of the galaxy samples consid- 



8 Note that, as a consequence, our notation differs from the one in 
Mandelbaum et al. (2006), Hirata et al. (2007), and Mandelbaum et al. 
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ered, the galaxy pairs used to compute ^ ot (r p , ft, z m ) consist of 
one galaxy out of the density sample tracing the galaxy distribu- 
tion and one galaxy out of the subsample with shape information 
used to trace the intrinsic shear. 

In Appendix A we derive an approximate procedure to 
model f? ot (r p , ft, Zm) to good accuracy. As we show there, 
, fl, Zm) can be obtained via a simple coordinate trans- 



-phot,- 



formation, see (A. 10), from the angular correlation function 



~2~7 



(4) 



where the angular gl power spectrum C g i is given in terms of the 
underlying three-dimensional power spectrum P s \ by the Limber 
equation 



Cgx(t;z\,Z2) 



X' 2 



(5) 



Here, x denotes comoving distance, integrated up to the comov- 
ing horizon distance ^hor- We have introduced the probability 
distributions of comoving distances p„ for the galaxy density 
sample and p e for the galaxy shape sample. They are related via 
Px(x\x(zd) = Px(z\z,i)dz/dx to the probability of a redshift z 
given the photometric redshift estimate I, . The latter can be ex- 
tracted from the histogram in Fig. 2 by a vertical section at zi- 
Note that we use these distributions, extracted from the 2SLAQ 
photometric versus spectroscopic relation and linearly interpo- 
lated, for both the shape-selected and the full number density 
sample as we find that their redshift distributions agree to good 
accuracy. This agreement is not obvious because the images 
of the galaxies selected for shape measurement need to have a 
certain minimum angular size. We trace the agreement of the 
two redshift distributions for the MegaZ-LRG sample back to a 
rough balance between the effect that galaxies at higher redshift, 
i.e. at larger distance, appear smaller due to the larger angular 
diameter distance, and the counteracting effect that for a given 
range of apparent magnitudes, galaxies at higher redshift are on 
average intrinsically brighter and thus have larger intrinsic phys- 
ical sizes, see also Sect. 2.5. 

Throughout, we will assume that the galaxy bias b g is scale- 
independent. Then P g i can be related to the matter-intrinsic 
power spectrum via P gl (k, z) = b s P S i (k, z), where we calculate 
Psi according to the linear alignment model (Catelan et al. 2001; 
Hirata&Seljak2004) 



P sl (k,z) = -A Ci p c 



(6) 



where D(z) denotes the linear growth factor, normalised to unity 
today. Equation (6) differs from the result derived by Hirata & 
Seljak (2004) in the dependence on redshift, with the latter ex- 
pression featuring an additional term (1 + z) 2 . In Appendix B 
we derive the correct scaling with redshift shown in (6); see also 
Hirata & Seljak (2010). We absorb the dimensions of the nor- 
malisation into the constant C\ which is set to C\ p a - ~ 0.0134, 
following Hirata & Seljak (2004) and Bridle & King (2007) 
who matched the amplitude of the linear alignment model to 
SuperCOSMOS observations at low redshift. This choice is 



(2010) in that these works used the term 'GI' for both the measured 
galaxy number density-ellipticity cross-correlations (gl in this paper) 
and the derived GI intrinsic alignment term. 
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Fig. 4. Three-dimensional correlation function models of a sample with 
the MegaZ-LRG photometric redshift error as a function of comov- 
ing line-of-sight separation IT and comoving transverse separation r p 
at Zm ~ 0.5. The galaxy bias has been set to 1.9 in all panels. Top panel: 
Galaxy clustering correlation (gg). Contours are logarithmically spaced 
between 1 (yellow shading) and 10~ 5 (violet shading) with three lines 
per decade. Centre panel: Galaxy number density-intrinsic shear corre- 
lations (gl). Contours are logarithmically spaced between 10~ 3 (yellow 
shading) and 10~ 6 (violet shading) with three lines per decade. Bottom 
panel: Galaxy-galaxy lensing (gG). For ease of direct comparison, the 
contours are encoded exactly like in the centre panel. Note that the 
galaxy-galaxy lensing signal is not symmetric around II = 0, in con- 
trast to the gg and gl terms. Also, it is negative, so that the modulus is 
plotted. For an illustration of the effect of photometric redshift scatter 
see Fig. A. 1. 



common but of no particular relevance for our study, and the 
normalisation is in principle arbitrary. Thus, we introduce a di- 
mensionless amplitude parameter A which is free to vary. If not 
stated otherwise, we use A — 1 to demonstrate our model in the 
following subsections. 
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The original derivation of the intrinsic alignment model re- 
quires the linear matter power spectrum in (6), see e.g. Appendix 
B, but following Bridle & King (2007) we use the full P s with 
nonlinear corrections, which provides satisfactory fits to existing 
data. Thus we refer to (6) as the nonlinear version of the linear 
alignment (NLA) model henceforth. The matter power spectrum 
is computed for our default cosmological model, a spatially flat 
ACDM universe with parameters £2 m = 0.25, cr 8 = 0.8, h = 0.7, 
and n s — 1 .0. The transfer function is calculated according to 
Eisenstein & Hu (1998) using £2b = 0.05, and nonlinear correc- 
tions are included following Smith et al. (2003). 

For illustration, in Fig. 4, centre panel, the predicted gl cor- 
relation function for the MegaZ-LRG sample (assuming align- 
ments consistent with those in SuperCOSMOS, and including 
photometric redshift errors) is plotted, for z m ~ 0.5 and assum- 
ing b g = 1 .9 which is roughly in agreement with the results ob- 
tained for the full MegaZ-LRG sample by Blake et al. (2007) 
and close to the value we determine, see Sect. 5.2. The signal 
is strongest around II = 0, but extends far out along the line- 
of-sight direction due to the photometric redshift scatter. The 
correlations have a maximum at r p ~ 0.5 /r'Mpc and decrease 
for larger r p due to the diminishing physical interaction between 
galaxies at large separation, and for small r p since the separation 
vector between pairs of galaxies gets close to the line-of-sight 
direction, see also Fig. A.l. 

Note that throughout this work we do not include redshift- 
space distortions in our modelling. Since for both spectro- 
scopic and photometric data we integrate the correlation func- 
tions over the line-of-sight separation out to at least 60/z'Mpc 
and 90/T 1 Mpc, respectively, redshift-space distortions should 
have a negligible influence on the integrated signals (see also 
the discussions in Padmanabhan et al. 2007; Blake et al. 2007; 
Mandelbaum et al. 2010). In addition, in the latter case, the red- 
shift space distortions should be subdominant compared to the 
size of the photometric redshift errors. 

4.2. Contribution of other signals 

Due to the photometric redshift scatter, contributions to galaxy 
position-shape correlations other than the gl term may become 
important. In the weak lensing limit, the measured ellipticity of a 
galaxy image is the sum of the intrinsic ellipticity and the grav- 
itational shear, while the number density is determined by an 
intrinsic term (whose two-point correlation is the usual galaxy 
clustering) plus modifications by lensing magnification effects. 
Hence, in terms of angular power spectra one can write (for de- 
tails see Bernstein 2009; Joachimi & Bridle 2010) 

Cne^zuzi) = Cgi(£;zi,a) + Cgo(*;zi,a) (7) 
+C mG (£;zuZ2) + C m i{£;z\,zi) , 

for each set of galaxy samples that is correlated. Apart 
from the gl signal, contributions from galaxy-galaxy lensing 
(gG), magnification-shear correlations (mG), and magnification- 
intrinsic correlations (ml) occur. 

If Z\ ~ Zi, the gl term is expected to dominate 9 , whereas 
galaxy-galaxy lensing is the dominant term in C ne {(; z\ , Zi) if 



9 It is instructive to bear in mind that the gl term, which we use to 
constrain the power spectrum P 6l , is restricted to pairs of galaxies at 
similar redshifts while the Gl signal, which is generated by P 6h is not 
limited to close pairs of galaxies and is actually strongest if galaxies 
are widely separated in redshift. The reason for these opposing redshift 
dependencies are the different weightings in the Limber equations, see 
(5) and (22) below. 



a number density and a shape sample at significantly different 
Z\ < Z2 are correlated. In addition, correlations between lensing 
magnification and gravitational shear can have a contribution, 
e.g. if a matter overdensity causes both tangential shear align- 
ment and an apparent boost in the number density of background 
galaxies. Likewise this overdensity could tidally align surround- 
ing galaxies and thus create correlations between magnification 
and the intrinsic galaxy shapes. 

All these additional signals are related to the three- 
dimensional correlation function via relations of the form (4), 
so that, to assess the importance of their contributions, it 
is sufficient to compare the angular power spectra. We re- 
strict the consideration to power spectra that correlate galaxy 
shapes and number densities at the same redshift (redshift auto- 
correlations), and that hence are representative of correlation 
functions at small line-of-sight separations. Note that for larger 
values of |I1|, the scaling of the different signals with redshift be- 
comes important, so that e.g. the amplitudes of the gG and mG 
signals relative to the gl term increase if bigger values for II max 
are chosen. 

The corresponding Limber equations of the additional sig- 
nals read (e.g. Joachimi & Bridle 2010) 



Cgotf;zi,z 2 ) = h A X ' 
Jo 



(8) 



X 



Cmo(f;zi,z 2 ) = 2(a-l) d x ' 
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where a — -d [log A^(> L)] /d [log L] is the logarithmic slope of 
the cumulative galaxy luminosity function of the density sample 
(e.g. Narayan 1989; Bartelmann & Schneider 2001). Moreover 
we have defined the lensing weight function 



qxtx, xi) 



3#oQm x 

2 c 2 aix) 



r 

•Jx 



dx' pAx'lXi) 



x -x 

X' 



(9) 



for x = {«, e), where a denotes the cosmic scale factor. 

To determine a, we calculate the cumulative distribu- 
tion of i band de Vaucouleurs magnitudes used for selec- 
tion of the MegaZ-LRG samples, and fit the slopes s = 
d [log iV(< z'dev)] /d/dev at the faint end of each distribution. The 
slope of the cumulative galaxy luminosity function is then given 
by a — 2.5 s (e.g. van Waerbeke 2010). We find a = 2.26 for the 
full MegaZ-LRG sample, a — 1 .29 for the low-redshift sample, 
and a = 3.04 for the high-redshift sample. 

In this comparison, we assume that the intrinsic alignment 
signal follows the corrected NLA model with the normalisa- 
tion from SuperCOSMOS. Since the strength of the intrinsic 
alignment signals for LRGs are expected to be significantly 
higher, see e.g. the results for SDSS LRGs by Hirata et al. 
(2007), this should be a conservative assumption. Note that the 
SuperCOSMOS normalisation employed in foregoing work was 
based on the redshift scaling of the II term given in Hirata & 
Seljak (2004). Since SuperCOSMOS has a mean redshift of 0.1, 
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Fig. 5. Top panel: Modulus of the angular power spectra of the differ- 
ent signals contributing to galaxy position-shape correlations. Number 
density-intrinsic correlations (gl) are shown in black, galaxy-galaxy 
lensing (gG) in red, magnification-shear correlations (mG) in green, 
and magnification-intrinsic correlations (ml) in blue. Solid curves cor- 
respond to the redshift auto-correlations at z = 0.5, i.e. close to the mean 
redshift of the low-redshift MegaZ-LRG sample, and dotted curves to 
the mean redshift of the high-redshift sample at z = 0.59. Bottom panel: 
Ratio of the aforementioned signals over the gl correlations, with the 
same coding of the curves as above. The grey region covers angular 
scales that do not contribute significantly to the correlation functions. 
Galaxy-galaxy lensing and possibly mG correlations yield a relevant 
contribution to number density-shape correlations besides the gl signal. 



the amplitudes of the II signals in the original and corrected ver- 
sion of the linear alignment model should differ by about a fac- 
tor of 1.1 (see also Appendix B). Thus we retain the normal- 
isation to SuperCOSMOS for the corrected version by setting 
A = 1.1 2 = 1.21. In addition, we choose a galaxy bias b % = 1.9, 
where the latter value is close to the actual fit results for MegaZ- 
LRG, see Sect. 5.2. The resulting angular power spectra for all 
four contributions to (7) are shown in Fig. 5. We plot the red- 
shift auto-correlation power spectra, using z\ — Zi = 0.5 and 
Z\ — Zi = 0.59 corresponding approximately to the mean red- 
shifts of the two MegaZ-LRG redshift-binned samples. 

For the photometric redshift accuracy of the MegaZ-LRG 
sample, the gl signal still clearly dominates the position-shape 
correlations. It has a slightly lower amplitude at z — 0.59 than 
at I = 0.5 due to the broader redshift distribution at the higher 
photometric redshift. To verify that it is indeed the width of the 
distribution, and not the shift of its mean redshift, that causes 
the depletion, we shift the redshift distribution at z = 0.4525 
to a mean of 0.6 and re-compute the gl signal which then has a 
similar, slightly higher amplitude compared to the gl correlations 
at I = 0.4525. The other signals are less affected by the width 
of the contributing redshift distributions since they depend on 
lensing and thus have a much broader kernel in the line-of-sight 
integration, see (8). 



The ml signal never attains more than a few per mil of the gl 
term and is hence irrelevant for our purposes. Due to the steep 
slopes of the galaxy luminosity functions, magnification-shear 
correlations can contribute up to 20 per cent of the gl term over 
a wide range of angular frequencies, and even become the dom- 
inant contamination of the gl signal at small I. However, due to 
the Bessel function Ji in the kernel of (4), contributions from 
small £ are largely suppressed in ^°\f p , H, z m ). The first maxi- 
mum of J2 is at (9 ~ 3, and the maximum angle probed by our 
analysis is 3.1 deg corresponding to r p - 60/i~'Mpc at z = 0.4, 
so that only angular frequencies I > 60 are important, as indi- 
cated by the grey region in Fig. 5. Still, the mG term could add 
of the order 10 % to the total signal under pessimistic assump- 
tions, which is in the regime of the expected model parameter 
errors, so that we include the mG signal into our modelling. We 
note that this contribution is likely to be even more relevant for 
future intrinsic alignment analyses of surveys with higher statis- 
tical power and/or less accurate photometric redshifts, provided 
that the faint-end slope of the luminosity function has similar 
steepness. 

Since we assume the NLA model, for which Pgi has the same 
angular dependence as the matter power spectrum, the differ- 
ences between the signals shown in Fig. 5 can only be due to the 
weights in the Limber equations (8). The gl, ml, and gG signals 
all have at least one term p(x) in the kernel which is thus rela- 
tively compact. Only the weight of the mG correlations features 
a product of lensing efficiencies (9) which smears out the infor- 
mation over comoving distance due to its width and in addition 
shifts the features in the projected power spectrum because q(x) 
peaks at half the source distance, see (9). Therefore the mG sig- 
nal has a different angular dependence, causing the waviness in 
the ratio mG/gl. 

Galaxy-galaxy lensing has a scale dependence that is similar 
to the gl term, thereby yielding a nearly constant contribution of 
about 30 %. Therefore we need to incorporate the gG term into 
our model, mainly affecting the amplitude of the model correla- 
tion function due to the almost constant ratio gG/gl. Note that 
contrary to the usual approach to galaxy-galaxy lensing studies, 
we have defined the correlation function such that radial align- 
ment produces a positive signal. Hence, the inclusion of the gG 
term into the model will increase the measured gl amplitude. 
The modulus of the three-dimensional gG correlation function 
is shown in Fig. 4, bottom panel. Due to the lensing contribu- 
tion, the gG correlation is not symmetric with respect to n = 0, 
even if the redshift distributions of the galaxy shape and density 
samples are identical. We note in passing that if the measure- 
ments included correlations between galaxies at largely different 
redshifts to which galaxy-galaxy lensing would be the dominant 
contribution, it would be possible to simultaneously measure the 
gG and gl signals, see e.g. Joachimi & Bridle (2010). However, 
such joint analyses are beyond the scope of the present work. 

The above statements hold only if the amplitude of the intrin- 
sic alignment signal is of the order found in the SuperCOSMOS 
survey, because we have used A ~ 1 in (6). If the contribution by 
intrinsic alignments were weaker, the importance of the gG and 
mG signals would further increase. However, Hirata et al. (2007) 
have demonstrated that LRGs show a strong intrinsic alignment 
signal at z ~ 0.3, so unless we find a strong decline of intrin- 
sic alignments with redshift, A ~ 1 should be a conservative 
assumption. 

We also consider galaxy clustering (hereafter gg) which will 
be used to determine the galaxy bias of the different samples. 
Since in the case of MegaZ-LRG the gg signal is affected in the 
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same way by the photometric redshift scatter as number density- 
shape cross-correlations, we proceed in exact analogy and com- 



pute the three-dimensional correlation function 
from (e.g. Hu & Jain 2004) 



phot 
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again by means of (A.l 1), where the angular power spectrum is 
related to the matter power spectrum via 
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Note that one of the redshift probability distributions is deter- 
mined from the shape sample because we use cross-correlations 
between the full and shape samples to measure galaxy clustering 
in this analysis, see Sect. 3. We show the three-dimensional cor- 
relation function of galaxy clustering in the top panel of Fig. 4. 
The strong spread of the gg signal along the line of sight demon- 
strates that in the case of the MegaZ-LRG sample, photometric 
redshift scatter and the corresponding effect of a truncation at 
large II when computing the projected correlation function has 
to be modelled with similar care as for the gl term. Since galaxy 
clustering produces a strong signal, we can safely neglect poten- 
tial contributions by lensing magnification effects in this case. 



4.3. Projection along the line of sight 

As in the spectroscopic case, the quantity that is actually com- 
pared to the data is the projected galaxy position-shape correla- 
tion function w g+ , obtained by integrating the three-dimensional 

correlation function ^ ot (f p ,fl, z m ) over Yl. In addition we take 
the average over a range of photometric redshifts z m which e.g. 
corresponds to the two redshift bins defined for the MegaZ-LRG 
sample, resulting in 

= f dtl f nV(z m ) e*°\r p ,n,z m ) , (12) 
J -n max J 

where the truncation at Il max , taken to be the same as for the data, 
has now been written explicitly. 

The average over z m contains the weighting 1V(z) over red- 
shifts as derived by Mandelbaum et al. (2010), which is given 
by 



nv(z) 



P 2 (z) 



X 2 (z)x'(z) 



dz 



P 2 (z) 



X 2 (z)x'(z) 



(12). A flux-limited sample like MegaZ-LRG misses faint galax- 
ies at high redshifts compared to a volume-complete sample, so 
these redshifts are down weighted accordingly by (13) in the av- 
eraging process. Note that this procedure of computing the three- 
dimensional correlation function for each redshift and subse- 
quently averaging over redshift with the weighting (13) is equiv- 
alent to our treatment of the data, where the correlation function 
was computed as a function of r p and II for all redshifts, thereby 
obtaining an average over the full range of redshifts covered by 
the sample. 

In the case of galaxy clustering, the projected correlation 
function w gg (r p ) is determined in exact analogy to (12). In prin- 
ciple the integral over ft in (12) should extend to infinity, see 
(A. 6), but it can be truncated at some maximum value fl max be- 
cause correlations between pairs of galaxies do not extend to 
infinite separation. However, the photometric redshift scatter for 
the MegaZ-LRG sample causes significant contributions to the 
line-of-sight integral in (12) from II > 300/i _1 Mpc, although 
signal-to-noise requirements enforce relatively small values of 

IT-max- 

For a consistency check on the modelling of the line-of-sight 
truncation of the correlation functions in case of the MegaZ- 
LRG sample, we compute the model and the observed corre- 
lation functions for fl max = 90/T 1 Mpc and n max = 180/z~'Mpc. 
For the number density-shape correlations as well as the galaxy 
clustering signal, we compare the ratios of the correlation func- 
tions with these two cut-offs in Sect. 5.1, finding good agreement 
between model and observational data. We use cut-offs in the 
signal integration along the line-of-sight at either 180/r'Mpc 
or 90/r'Mpc for the fits to the full MegaZ-LRG sample and 
find consistent results for the intrinsic alignment fit parame- 
ters with errors of the same order, see Sect. 5.3. The signals for 
w gg and w g+ both have similar signal-to-noise when truncating 
at these two values of n max . The correlation functions for the 
two MegaZ-LRG redshift-binned samples have been cut off at 
n max = 90/r'Mpc throughout. 

In addition to the photometric MegaZ-LRG sample, we will 
also reconsider spectroscopic samples from SDSS. As discussed 
before, the line-of-sight truncation can be ignored in the case of 
spectroscopic data, so that the projected correlation function is 
simply given by 

dzOV(z) J o J 2 (k±r p ) P a (k ± ,z) ,(14) 

where (A. 6) and the same redshift averaging procedure as in (12) 
were used. Similarly, one obtains for the spectroscopic galaxy 
clustering signal (e.g. Hirata et al. 2007) 

//"*00 17 7 
dzOV(z) J o -^-±J (k ± r p )P s (k ± ,z). (15) 



where p(z) is in this case the unconditional probability distri- 
bution of photometric redshifts determined from the MegaZ- 
LRG sample 10 (or its redshift-binned subsamples). Here, x'(z) 
denotes the derivative of comoving distance with respect to red- 
shift. Note that the denominator x 2 (z)x'(z) in (13) is propor- 
tional to the derivative of the comoving volume V com with re- 
spect to redshift. Equation (13) can be illustrated by considering 
a volume-limited sample for which p(z) = dV com /dz holds. Then 
'W(z) = p(z), as expected for a simple average over redshift in 



10 Note that using the photometric redshift distribution obtained from 
the 2SLAQ verification sample instead has no significant influence on 
the models. 



4.4. Fitting method 

We perform the fits to the data via weighted least squares min- 
imisation, using the reduced x 2 at the minimum to quantify the 
goodness of fit. The data are compared to signals computed for 
the NLA model (6) or a more flexible variant (19) introduced 
in Sect. 5. To obtain confidence regions, we compute the likeli- 
hood according to L oc exp(-^ 2 /2), i.e. assuming Gaussianity. 
The posterior probability in parameter space is computed via 
Bayes' theorem, using a top-hat prior that is truncated outside 
the regime where the likelihood deviates substantially from zero, 
see Sect. 5.5 for details. When doing three-parameter fits, we 
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Fig. 6. Correlation coefficients r con between different transverse separation bins of the correlation functions w u (left panel) and w g+ (right panel) 
for the redshift-binned MegaZ-LRG samples. Shown is the modulus of r corr for the 10 r p bins correlated with the smallest bin considered (r p = 
0.44 Mpc, solid lines), the smallest bin used for the fits (r p = 7.67 /r'Mpc, dot-dashed lines), the centre bin used for the fits (r p = 19.25 /r'Mpc, 
dashed lines), and the largest bin (r„ = 48.30 /r'Mpc, dotted lines). Black curves correspond to the lower redshift bin, red curves to the higher 



redshift bin. Note that the black and red lines for the correlations of w„. 



with r p = 0.44 h~ 



Mpc nearly coincide. 



marginalise in each case over the hidden parameters, lcr and 2cr 
confidence contours are then defined by the regions containing 
68.3 % or 95.4 % of the (marginalised) posterior. 

Since we have noisy jackknife covariances obtained from a 
finite number of realisations, the inverse of these covariances, 
required for the % 2 , is biased (Hirata et al. 2004; Hartlap et al. 
2007). We employ the corrected estimator for the inverse covari- 
ance presented in Hartlap et al. (2007), 



(Cov)- 1 = 



n — d — 2 
n- 1 



(Cov)"' = T (Coy) 



(16) 



where d is the dimension of the data vector and n the number of 
realisations used to estimate the covariance. Equation (16) was 
derived under the assumption of statistically independent data 
vectors with Gaussian errors. For the SDSS samples (d = 10, n = 
50) we find T ~ 0.776, and for MegaZ-LRG {d = 10, n = 256) 
T ~ 0.957, the latter result being in excellent agreement with 
the results obtained from the simulations described in Appendix 
D of Hirata et al. (2004). 

To study the characteristics of the covariances, we compute 
the correlation coefficient between different transverse separa- 
tion bins, 



^"corr(/ 'p,a-> *~p,b) ~ 



Cov(r Pta ,r Ptb ) 



- v /Cov(r Pjfl , r PtCI ) Cov(r Ptb , r p , h ) 



(17) 



In Fig. 6 we have plotted r corr of both w gg and w g+ for the two 
MegaZ-LRG samples at low and high redshift. While w g+ decor- 
relates quickly with only moderate correlation between neigh- 
bouring bins on the largest scales, w gg features strong positive, 
long-range correlation particularly on the larger scales that are 
used for the fits. The correlation coefficients have similar val- 
ues for the two redshift bins, with the z < 0.529 bin showing 
slightly higher correlation in most cases. For the SDSS samples 
we find a similar correlation structure for w g+ , but much weaker 
correlations in w„ s . 

BB 

The difference in correlation length between w gg and w g+ is 
caused by the different kernels in the Hankel transformation be- 
tween power spectrum and correlation function, see (4) and (10). 



Since Jiix) decreases faster than Jq(x) for increasing x, we ex- 
pect w gg to generally feature stronger correlations. A given trans- 
verse separation r p between galaxies is observed under a smaller 
angle if these galaxy pairs are located at higher redshift, and it 
is this angle which enters the argument of the Bessel functions. 
Therefore the correlation present in w gg is more pronounced in 
the MegaZ-LRG samples, which are at considerably higher red- 
shift than the other SDSS samples. 



5. Results 

5.1. Scaling with line-of -sight truncation 

To test whether the data and the model show the same behaviour 
when varying n max , we compute both w g+ and w gg for the 
MegaZ-LRG sample according to (12) for n max = 90/i~'Mpc 
and n max = 180/; _1 Mpc. Then we compare the ratio of w g+ 
with cut-off n max = 90/z~'Mpc over w g+ with cut-off n max = 
180/; _1 Mpc (and likewise for w gg ) obtained from the model to 
the corresponding ratio computed from the observations. Since 
the projected correlation functions with the different cut-offs are 
strongly correlated, we compute the errors on the ratio again via 
jackknifing. Note that due to these correlations the actual errors 
on the ratio are significantly smaller than if one assumed them 
to be independent. Note furthermore that the ratio also inherits 
a significant correlation between different transverse separations 
from the individual projected correlation functions, in particular 
for w gg , see Fig. 6. 

In Fig. 7 we have plotted the ratios of the projected corre- 
lation functions with the different cut-offs. The model predic- 
tion for the ratio of the galaxy clustering signals is in fair agree- 
ment with the data, yielding a loss of about 40 % when reduc- 
ing LI max from 180 /r'Mpc to 90 /r'Mpc. On scales where the 
galaxy bias becomes nonlinear, indicated by the grey shaded re- 
gion, bias effects may not cancel in the ratio anymore, so that 
deviations from the model prediction are expected. On scales 
r p > 20 /r' Mpc one observes an apparently significant trend in 
the data to fall below the model prediction. However, note that 
w gg at large transverse separations features very strong cross- 
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Fig. 7. Effect of the cut-off in IT in the projection of the three- 
dimensional correlation functions along the line of sight. Shown is 
the ratio of the projected correlation function computed for n max = 
90/r'Mpc over the correlation function obtained with n max = 
180/7~'Mpc, for both the galaxy clustering signal (gg, in black) and 
number density-shape correlations (g+, in red). Points are computed 
from the MegaZ-LRG data, using the full range in redshifts. Note that 
the black points have been slightly offset horizontally for clarity. The 
black line is obtained from the model for w gg . The red hatched region 
comprises the range of ratios for w g+ with different relative strengths of 
the galaxy-galaxy lensing contribution. The red solid line indicates the 
ratio resulting for the best-fit intrinsic alignment amplitude determined 
in Sect. 5.3. Note that the error bars at different transverse separations 
are correlated, in particular for w gg at large r p , see Fig. 6. 



correlations between the data points, and this property is inher- 
ited by the ratio. A fit of the model ratio to the observed ratio 
for the five data points with r p > 6/z~'Mpc yields a reduced^ 2 
of 1.75, corresponding to a /5-value of 0.12, and hence data and 
model are still marginally consistent. 

The prediction for the ratio of the position-shape correlation 
functions with cut-offs at 90/i~'Mpc and 180/z~'Mpc, respec- 
tively, depends on the relative strength of the galaxy-galaxy lens- 
ing contribution with respect to the gl signal (we neglect the mG 
contribution whose effect should be very small in this case). We 
shade the possible range of ratios in Fig. 7 between the lowest 
ratio resulting for a negligible gG term which almost coincides 
with the ratio for galaxy clustering, and the highest ratio result- 
ing for the set of parameters we used in Sect. 4.2 which we take 
as a conservative lower limit on the intrinsic alignment signal. 
We also show the curve obtained with the best-fit intrinsic align- 
ment amplitude from the fits in Sect. 5.3 below, which is in good 
agreement with the data, yielding a ratio of 0.64. In this case 
the least-squares analysis results in a reduced x 2 of 0.77 and a 
p-value of 0.57. 

Note furthermore that both models and data are consistent 
with the fact that the loss of signal due to the smaller cut-off 
in Tl is roughly constant in transverse separation. The general 
agreement of the observed and modelled ratios confirms that we 
model the effect of n max on photometric redshift data correctly. 
Besides, it supports our use of the 2SLAQ photometric redshift 
error distribution despite the slightly different apparent magni- 
tude limits of the sample, as discussed in Sect. 2.4. 



5.2. Galaxy bias 

To relate the observed galaxy number density-intrinsic corre- 
lations (plus the corrections due to galaxy-galaxy lensing and 
magnification-shear correlations) to the matter-intrinsic correla- 
tions that generate intrinsic alignments, the galaxy bias b g for 
the density tracer sample needs to be measured. As described 
in Sect. 3, we compute a galaxy clustering signal that represents 
the cross-correlation between the density tracer sample and the 
sample used to trace the intrinsic shear. Given that the latter is 
a subset of the former with nearly the same properties (redshift 
and luminosity distributions, see Fig. 3), we assume that the two 
have the same galaxy bias. Thus, we compute b g from this galaxy 
clustering signal, assuming a linear bias model, but using the full 
matter power spectrum which should extend the validity of the 
fits into the quasi-linear regime (see also Hirata et al. 2007, who 
test several methods to determine the galaxy bias in a similar 
context). Note that all our considerations rely on the hypothesis 
that we have assumed the correct cosmological model, in partic- 
ular cr 8 = 0.8. 

The redshift averaging and the projection along the line of 
sight of Wgg is performed according to (12) for the photomet- 
ric MegaZ-LRG samples and following (15) for the SDSS LRG 
samples. We do not repeat the bias measurement for the SDSS 
L3 and L4 samples but adopt the values determined by Hirata 
et al. (2007), rescaled to our value of cr 8 by employing b g oc <x~ 1 , 
which results in b g = 1 .04 and b g = 1 .01, respectively. Note that 
the assumption of that same bias, despite the use of different 
colour cuts for the intrinsic shear tracers, is acceptable because 
the bias we need is that of the density tracer sample, which has 
not changed. To all model projected correlation functions w gg , 
we add a constant C as a further fit parameter to account for the 
undetermined integral constraint on the numerator of the estima- 
tor (3) due to the unknown mean galaxy number density (Landy 
& Szalay 1993; see also Hirata et al. 2007), i.e. 



w^ e \r p , b g , C) = WggOv, b g ) + C, (18) 

where w gg is given by (15) or the analogue of (12) in the case 
of the MegaZ-LRG samples. Note that we have made the depen- 
dence on the galaxy bias explicit in the foregoing equation. 

For the fit we discard scales r p < 6 /z~'Mpc, i.e. the five data 
points at the smallest r p where the assumption of a linear bias 
is expected to break down (Tasitsiomi et al. 2004; Mandelbaum 
et al. 2006). For the MegaZ-LRG sample, there is one additional 
nuisance in this modelling, which is the 5 % stellar (M star) con- 
tamination fraction. As shown in Mandelbaum et al. (2008), the 
imposition of the apparent size cut that is needed for a robust 
galaxy shape measurement is sufficient to remove this contam- 
ination to within 1%. Thus, the galaxy clustering signal as de- 
fined (a cross-correlation between the shape and density sam- 
ples) is diminished by a single power of the contamination frac- 
tion /contam = 0.05. The bias determined by the fits is actually 
\/l ~ /contain b g , so we must correct it upwards to account for 
that. Then, since w g + is reduced by a factor of 1 - /contain, instead 
of dividing by j ust b g to get to w 5+ , we divide by ( 1 - / coltam )b„. 

In Fig. 8, we show the projected correlation functions w gg for 
the two MegaZ-LRG redshift bins and the two SDSS LRG red- 
shift bins. Note that the SDSS LRG samples have not been split 
further into luminosity bins because the full SDSS LRG sam- 
ple, divided into the two redshift bins, is used to trace the galaxy 
number density field. In each case, we also plot the best-fit mod- 
els, indicating that the model is a good description of the data on 
scales where nonlinear bias is not important. At smaller trans- 
verse separations (the grey region), which have been excluded 
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Table 2. Galaxy bias b„ for the different galaxy samples used. 



sample 


h 


2 

^red 


: — ; T\ — 

p(>x l ) 


MegaZ-LRG, all z, (90) 


1.94 ±0.03 


0.21 


0.89 


MegaZ-LRG, all z, (180) 


1.91 ±0.03 


0.37 


0.77 


MegaZ-LRG, z < 0.529, (90) 


1.87 ±0.04 


0.22 


0.88 


MegaZ-LRG, z > 0.529, (90) 


2.09 ± 0.04 


2.51 


0.06 


SDSS LRG, z < 0.27 


1.88 ±0.10 


0.82 


0.48 


SDSS LRG, z > 0.27 


1.89 ±0.07 


0.97 


0.41 


SDSS Main L4 red 


1.04 






SDSS Main L3 red 


1.01 







Notes. Shown is the lcr-error on b g , marginalised over the additive con- 
stant C, the reduced^ 2 for 3 degrees of freedom, and the corresponding 
p- value. The values for the SDSS L3 and L4 samples have been derived 
from the results by Hirata et al. (2007). All results assume <x 8 = 0.8. 
The numbers in parentheses indicate n max in units of /r'Mpc for the 
MegaZ-LRG samples. 



from the fits, the data have increasingly larger positive offsets 
with respect to the model, caused by nonlinear clustering effects. 

The best-fit values for b g , marginalised over C, are sum- 
marised in Table 2. We find good agreement within the Icr lim- 
its for the best-fit galaxy bias, determined for the different Il max 
in the MegaZ-LRG data, again confirming that we are correctly 
modelling the truncation in the line-of-sight projection. Splitting 
the MegaZ-LRG data into two redshift bins at z = 0.529, we 
obtain a stronger bias for the bin at higher redshift. This finding 
is expected, as the bin with z > 0.529 contains on average sig- 
nificantly more luminous galaxies that are more strongly biased, 
see Fig. 1. Only the high-redshift MegaZ-LRG sample yields a 
reduced x 1 that significantly exceeds unity which we trace back 
to the strong correlations between errors as the plot in Fig. 8 sug- 
gests a good fit. Indeed the reduced x 1 drops well below unity 
if we repeat the fit ignoring the off-diagonal elements in the co- 
variance. 

We compare our results with the galaxy bias obtained by 
Blake et al. (2007) who also studied MegaZ-LRG, albeit with 
slightly different selection criteria. They used the cuts /d e v ^ 
19.8 and d± > 0.55 throughout, as well as additional star- 
galaxy separation criteria that reduced the stellar contamination 
to 1.5 %. The different selection criteria hinder direct compari- 
son, e.g. the Blake et al. (2007) criteria (driven mostly by the /d e v 
cut) shift the r band absolute magnitude range 0.15 mag brighter. 
Hence, we expect that the galaxy bias obtained from our analysis 
should be smaller, and indeed, after rescaling the bias given in 
Table 2 of Blake et al. (2007) to cr 8 = 0.8, we find b g = 1.89 and 
b s = 2.10 for the two redshift slices roughly coinciding with our 
MegaZ-LRG low-redshift sample, and b g = 2.18 and b g = 2.44 
for the two redshift slices closer to our high-redshift sample. 

The SDSS LRG samples yield a similar galaxy bias com- 
pared to the full MegaZ-LRG sample, with no significant evolu- 
tion in redshift. Given that the SDSS LRG galaxies have on av- 
erage a higher luminosity and are located at considerably lower 
redshift, this finding hints at a stronger bias in the past for galax- 
ies at fixed luminosity. Using again the fact that the bias scales 
as b g oc err 1 , our findings for the SDSS LRG samples can be 
compared to the results for the equivalent bias model in Hirata 
et al. (2007) who use cr 8 = 0.751. Rescaling the values of Table 
2 to this value of erg, we get b g = 2.00 + 0.11 for the low-redshift 
sample and b g = 2.01 + 0.07 for the high-redshift sample. These 
values agree (within Icr) with b g = 2.01 ± 0.12 for z < 0.27 
and b g = 1.97 + 0.07 for z > 0.27 as found by Hirata et al. 
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Fig. 8. Projected correlation function w SB as a function of comoving 
transverse separation r p . Top panel: For the SDSS LRG sample with 
redshifts smaller than 0.27 (black) and with redshifts larger than 0.27 
(red). Bottom panel: For the MegaZ-LRG sample with photometric red- 
shifts smaller than 0.529 (black) and with photometric redshifts larger 
than 0.529 (red). Note that the red points have been slightly offset hor- 
izontally for clarity, and that the error bars are strongly correlated. In 
addition we show the best-fit models as black and red curves, respec- 
tively. Only the data points outside the grey region have been used for 
the fits to avoid the regime of nonlinear bias. 



(2007). Note that the latter analysis used a narrower range in 
transverse separation with r p = 7.5 - 47 /r'Mpc compared to 
r p = 6 - 60 /i~'Mpc considered in this work. 

5.3. Intrinsic alignment model fits to individual samples 

With the galaxy bias in hand, we can now proceed to fit mod- 
els of intrinsic alignments to w s+ . The NLA model features a 
single free parameter for the amplitude, A. Within the physical 
picture of this model, the amplitude quantifies how the shape of a 
galaxy responds to the presence of a tidal gravitational field. It is 
likely that this response depends on the galaxy population under 
consideration, and thus possibly features an additional evolution 
with time and hence redshift dependence (on top of the one in- 
herent to the NLA model), and a variation with galaxy luminos- 
ity. Therefore we will investigate a more flexible prescription for 
the gl power spectrum in Sect. 5.5. 

In this section we use (6) as the intrinsic alignment 
model, with the single fit parameter A. We keep the original 
SuperCOSMOS normalisation, i.e. C x p a ~ 0.0134. To allow 
for a comparison with foregoing work, we also present some 
fits with models based on the NLA version with the redshift de- 
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Fig. 9. Projected correlation function w g+ as a function of comoving 
transverse separation r p for different MegaZ-LRG subsamples. Top 
panel: Shown is w g+ for the MegaZ-LRG sample with the full range in 
redshifts. The black solid curve corresponds to the best-fit model when 
only varying the amplitude A, without contributions by galaxy-galaxy 
lensing. The dotted black line is obtained by using the linear matter 
power spectrum instead of the full power spectrum including non-linear 
corrections, and identical model parameters otherwise. Bottom panel: 
Same as above, but for the MegaZ-LRG sample split into two photo- 
metric redshift bins, where results for z < 0.529 are shown in black, and 
for z > 0.529 in red. Dotted lines again correspond to models computed 
from the linear matter power spectrum. The point for the z < 0.529 sub- 
sample at r p as 1 /i~'Mpc is negative and thus not shown. Note that the 
red points have been slightly offset horizontally for clarity, and that the 
error bars are correlated. Only the data points outside the grey region 
have been used for the fits. 



pendence given in Hirata & Seljak (2004). Note that all intrin- 
sic alignment models applied in this work have a fixed depen- 
dence on transverse scales. Since the assumption of a linear bias 
also enters the model, we again limit the parameter estimation 
to scales r p > 6/z~ 1 Mpc. Note that we do not explicitly prop- 
agate the errors on the galaxy bias determined in the foregoing 
section through to the uncertainty on intrinsic alignment param- 
eters, as they are marginal compared to the measurement error 
in w g+ (which is dominated by shape noise). 

In Fig. 9, the projected correlation functions for the full 
MegaZ-LRG sample as well as for the two MegaZ-LRG red- 
shift bins, split at z = 0.529, are plotted. The fit results for A 
are presented in Table 3. On the scales usable for the fit, the 
best-fit gl model, which is also plotted in Fig. 9 for each case, 
traces the data points well with reduced ^-values below one, 
whereas for r p < 1 /t 'Mpc points lie several cr above and below 
the model curve, possibly indicating strongly nonlinear effects. 



Fig. 10. Projected correlation function w„ + as a function of comoving 
transverse separation r,, for different SDSS samples. Top panel: Shown 
is w g+ for the red SDSS L4 (in black) and L3 (in red) samples. The 
curves correspond to the best-fit models when only varying the ampli- 
tude A. Bottom panel: Same as above, but for the SDSS LRG medium 
brightness sample split into two redshift bins, where results for z < 0.27 
are shown in black, and for z > 0.27 in red. The point for the z < 0.27 
subsample at the smallest r p is negative and thus not shown. Note that 
the red points have been slightly offset horizontally for clarity, and that 
the error bars are correlated. Only the data points outside the grey region 
have been used for the fits. 



The nature of these deviations is unknown, but since they occur 
on scales near the virial radius of LRGs, one may hypothesise 
that at these ranges of r p , complicated dependencies on the tidal 
field or a change in the intrinsic alignment mechanism play a 
role. Moreover we find very good agreement between the best-fit 
amplitudes obtained for the full MegaZ-LRG sample with differ- 
ent values of Il max . 

In addition, we show in Fig. 9 models for w„+ that have been 
calculated using the linear matter power spectrum instead of the 
nonlinear one in (6), holding all other model parameters fixed. 
As expected, the signals for linear and nonlinear power spectrum 
coincide on the largest scales, but already at r p ~ 10/i~'Mpc, 
w g+ computed from linear theory drops below the correlation 
function that includes nonlinear clustering and yields a worse 
fit to the data in case of the full and the high-redshift MegaZ- 
LRG sample. Thus, although our analysis is still restricted to 
relatively large scales, non-linear effects in the intrinsic align- 
ment of galaxies must be taken into account. 

We also perform the analysis on w g+ for the SDSS LRG data, 
which is divided into three luminosity bins in addition to the two 
redshift bins split at z = 0.27, see Table 1. As redshifts are de- 
termined spectroscopically in this case, we use (14) to compute 
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Table 3. lcr confidence limits on the amplitude A of the gl correlation function for each sample. 



sample g (L)/L A p(> X z ) A (HS04) 

MegaZ-LRG, all z, (90) 0.54 0.96 4 - 52 -°a 6 64 °- 04 1 00 1 - 98 -028 

MegaZ-LRG, all z, (180) 0.54 0.96 4.51™* 0.36 0.84 1.98 

i +o:86 



0.28 



MegaZ-LRG, z < 0.529, (90) 0.49 0.87 5.3L"™ 0.47 0.76 2AQ -nll 

MegaZ-LRG, z > 0.529, (90) 0.59 1.05 3.53^ 0.58 0.68 1.46™* 



MegaZ-LRG, all z, g - i cut 0.53 0.99 5.23 + " : " 0.26 0.90 2.29 + °^ 

MegaZ-LRG, z < 0.529, g - i cut 0.49 0.91 5.91+™ 0.17 0.95 2.67 + ° : « 

MegaZ-LRG, z > 0.529, g - i cut 0.59 1.08 4.53 + \f 6 0.20 0.94 1.87 + g : |f 

SDSS LRG, z < 0.27, faint 0.21 1.06 4.51+j™ 1.30 0.27 3.03+jlj 

SDSS LRG, z > 0.27, faint 0.32 1.07 7.67!™ 1.82 0.12 4.44+j : ^ 

SDSS LRG, z < 0.27, med 0.22 1.50 9.98+j : | 0.51 0.73 6.77+| : °? 

SDSS LRG, z > 0.27, med 0.31 1.50 10.03+j*? 2.03 0.09 5.84^ : | 

SDSS LRG, z < 0.27, bright 0.22 2.13 12.931™ 1.15 0.33 8.77 + ™ 

SDSS LRG, z > 0.27, bright 0.31 2.12 16.09+!™ 2.07 0.08 9.35 + f 

SDSS L4 red 0.10 0.33 1.20 + ^g 0.80 0.52 L00 + * : | 

SDSS L3 red 0.07 0.14 3.61 + ^ 0.83 0.51 3.16 + ] g 



sample 


A (w/o gG+mG) 


^red 


p(>x z ) 


A (w/o mG) 


MegaZ-LRG, all z, (90) 
MegaZ-LRG, all z, (180) 
MegaZ-LRG, z < 0.529, (90) 
MegaZ-LRG, z > 0.529, (90) 


A 1 1 +0.64 
^- 11 -0.64 

3.95+0.63 
4 (j2+8$> 
3 06™^ 


0.03 
0.36 
0.47 
0.58 


1.00 
0.84 
0.76 
0.68 


4 47+0.64 

4 45+0.64 

5 26™^ 
o 47+0.57 



Notes. Top section: Fits results for all samples considered in this work. For reference, we give the means (z) and (L) /Lq for each sample in the 
second and third column. The fourth to sixth columns contain the best-fit intrinsic alignment amplitude, and the corresponding reduced x 1 an d 
p- values of the fit for 4 degrees of freedom. To facilitate comparisons with preceding work, we also fit the amplitude of the version of the NLA 
model based on Hirata & Seljak (2004), shown in the last column. Note that these fits produce the same^ 2 as the foregoing ones. We also show 
the fit results for MegaZ-LRG samples with a colour cut as introduced in Sect. 5.4. The numbers in parentheses indicate II max in units of /r'Mpc 
for the MegaZ-LRG samples. Bottom section: Amplitude fits to MegaZ-LRG samples neglecting the contributions by galaxy-galaxy lensing (gG) 
and magnification-shear cross-correlations (mG). The rightmost column contains best-fit amplitudes when including the gG but not the mG signal. 
Due to the very similar r,,-dependence of the gl, gG, and mG signals, the ^-values of these fits are almost identical. 



w g +. For reasons of simplicity we use the redshift distribution 
for all SDSS LRG luminosities combined because we find that 
employing the individual distributions for the faint, medium, and 
bright subsamples instead leads only to sub-per cent changes in 
the signal on all scales. The resulting correlation functions and 
their best-fit models are shown in Fig. 10, and the resulting pa- 
rameter constraints on A listed in Table 3. 

In Table 3 we additionally present constraints on A using the 
NLA model with the redshift dependence as derived by Hirata 
& Seljak (2004). If the redshift distribution of a galaxy sample 
is sufficiently compact, the amplitudes of the two NLA mod- 
els considered are approximately related by a factor (1 + (z)) 2 , 
caused by the different redshift dependencies. Our findings for 
the SDSS LRG samples are compatible with the results of the 
power-law fits by Hirata et al. (2007), yielding a maximum in- 
trinsic alignment amplitude of close to ten times that found in 
SuperCOSMOS, for the bright high-redshift SDSS LRG sam- 
ple. In contrast, the resulting values for A using the SDSS L3 
and L4 samples are small for both NLA models, the constraints 
being consistent with zero at the 2cr-level. 

By default we include both galaxy-galaxy lensing and 
magnification-shear correlations in the modelling for the photo- 
metric redshift MegaZ-LRG data, but the bottom panel of Table 
3 also lists results for A when dropping either the mG term only 
or both additional contributions. Since the gG and mG signals 
yield a negative contribution to w g+ , a lower amplitude A than 
in the case including these contributions is needed to get a good 



fit to the data. Dropping the mG term causes a drop in A for all 
samples which is below the 2 % level and hence much smaller 
than the lcr error on the amplitude. Since the intrinsic alignment 
amplitude is about three to four times higher than assumed in the 
prediction of Sect. 4.2, which yielded a maximum mG contribu- 
tion of the order 10 % on relevant scales, effects at the per-cent 
level by the mG signal are indeed expected. 

The change in amplitude when discarding all additional sig- 
nals ranges between 7 % for the low-redshift MegaZ-LRG sam- 
ple and 13 % for the high-redshift sample. As Fig. 5 suggests, 
the relative contribution of the gG to the gl signal is approxi- 
mately constant as a function of redshift for MegaZ-LRG if the 
intrinsic alignment amplitude is fixed. Since we obtain a larger A 
for the low-redshift than for the high-redshift sample, the contri- 
bution by galaxy-galaxy lensing is correspondingly smaller for 
z < 0.529. In Sect. 4.2 we calculated a 30 % gG/gl ratio for 
A = 1.21, which is again in good agreement with our fit results. 

It is interesting to note that our default analysis yields nearly 
perfect agreement between the best-fit values for A obtained 
from the full MegaZ-LRG samples with cut-offs at 90/r'Mpc 
and 180/i~'Mpc, respectively, but that, when excluding galaxy- 
galaxy lensing, one observes a moderate discrepancy in the fit 
results. The galaxy-galaxy lensing increases for larger line-of- 
sight separation of the galaxy pairs correlated, whereas the gl 
term diminishes, so that the subsample with n max = 180 /r'Mpc 
is affected more strongly. This finding again confirms that we are 
modelling the effect of photometric redshift scatter and the trun- 
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Fig. 11. Constraints on the amplitude A of the intrinsic alignment 
model. The black solid curve corresponds to using the full MegaZ- 
LRG sample. The constraints from the individual MegaZ-LRG redshift 
bins are shown as red dashed lines (z > 0.529) and orange dotted lines 
(z < 0.529). For comparison we also show the constraints on A for the 
SDSS L4 sample as green solid line and for the SDSS LRG medium 
luminosity samples with z < 0.27 (dark blue dotted line) and z > 0.27 
(purple dashed line). 
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cation of signals at large II correctly. Besides, both additional 
signals that contribute to galaxy position-shape correlations have 
a dependence on transverse separation that is very similar to the 
one of the gl term. Thus only the amplitude A is affected while 
the goodness of fit remains almost unchanged when the gG and 
mG terms are included. 

In general, the dependence on r p given by the NLA model 
describes the data reasonably well, yielding reduced^ 2 values of 
order unity (or sometimes below). Only the high-redshift SDSS 
LRG samples tend to ax 2 that significantly exceeds unity, which 
is caused by an excess signal around r p = lO/^'Mpc of un- 
known origin, as can be seen in Fig. 10, bottom panel. The p- 
values remain above, but are close to, the significance level of 
0.05. 

In all cases except the two SDSS Main samples the in- 
trinsic alignment amplitude is higher than for the original 
SuperCOSMOS normalisation, which would correspond to A — 
1.21 in the corrected linear alignment model and A = 1 when 
fitting the original version. The SDSS Main samples are at sim- 
ilar redshift to SuperCOSMOS, and consistent with the results 
by Brown et al. (2002), although only at the 2cr level in the case 
of the L3 sample which prefers a higher amplitude (but note the 
possible excess signal at r p ~ lO/^'Mpc in the top panel of 
Fig. 10). 

As is obvious from the compact and mutually inconsis- 
tent posterior probabilities for A shown in Fig. 11, the differ- 
ent galaxy samples are inconsistent with an intrinsic alignment 
model that has only A as a free parameter. The SDSS LRGs 
span a very similar and relatively short range in redshifts, so 
that no strong evolution with redshift is expected in these sub- 
samples. Then it is evident from the fit results in Table 3 that the 
intrinsic alignment amplitude increases with galaxy luminosity, 
with the brightest sample attaining a high amplitude of A 16. 
Moreover, despite a mean luminosity that is 20 % higher than 
that of the low-redshift MegaZ-LRG sample, the high-redshift 
MegaZ-LRG sample has a smaller amplitude parameter A, in- 
dicative of a decrease of the intrinsic alignment amplitude with 
redshift beyond the redshift dependence inherent to the NLA 
model. However, note that the amplitudes of the two samples are 



Fig. 12. Redshift, colour, and magnitude properties of all samples used 
in this paper. Top panel: k + e-corrected (to z = 0) M r used to define 
the luminosities relative to L , versus redshift z. Galaxies from MegaZ- 
LRG are shown in green, those from the SDSS LRG samples in blue, 
and galaxies from the red SDSS Main L4 (L3) sample in red (black). 
To avoid confusion, a subset of the points from each sample was used. 
Bottom panel: The colour-redshift relation. The colour coding is the 
same as above. As shown, the full MegaZ-LRG sample is the only one 
with significant contributions below 2 (g - i) = 1.7. Thus, we define the 
cut sample (pink points) using all MegaZ-LRG galaxies that are redder 
than this value, for consistency with the other samples. 

still consistent with each other at about the lcr level, so that the 
NLA model is still in agreement with these findings for MegaZ- 
LRG. 

5.4. Compatibility of the different samples 

As has been shown in Hirata et al. (2007) and many other pa- 
pers, red and blue galaxies behave differently with respect to in- 
trinsic alignments. Thus, in order to obtain combined constraints 
from these samples and to understand their results in some uni- 
fied way, we now address the sample definitions as regards the 
resulting colour properties in greater detail. 

Wake et al. (2006) performed a comparison of the SDSS 
spectroscopic LRG and the MegaZ-LRG sample definition, with 
the goal of creating a subset of the MegaZ-LRG sample that 
would pass the SDSS LRG joint colour-magnitude cuts if it were 
shifted to redshift z = 0.2. The reason for this choice of compar- 
ison redshift is that, with the MegaZ-LRG sample concentrated 
at z = 0.55, this difference in redshift corresponds to shifting 
the SDSS filters over by one, and thus the ^-corrections are less 
prone to systematic error. As shown in their comparison, when 
using the g - i colour shifted to z = 0.2, only ~ 30 % of the 
MegaZ-LRG galaxies would pass the joint colour-magnitude cut 
of the SDSS spectroscopic LRGs. 

However, this cut is not what we want to impose on our sam- 
ple. The reason for this is that, as shown in Wake et al. (2006), 
Fig. 3, many of the MegaZ-LRG galaxies that are excluded are 
at the faint end, where the SDSS spectroscopic LRG cut starts to 
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Fig. 13. Projected correlation function w g+ as a function of comoving 
transverse separation r p for the MegaZ-LRG subsamples with the cut 
02 (g — i) > 1.7 imposed. Top panel: Shown is w i+ for the MegaZ-LRG 
sample with the full range in redshifts. The black solid curve corre- 
sponds to the best-fit model when only varying the amplitude A, without 
contributions by galaxy-galaxy lensing and magnification-shear cross- 
correlations. The dotted line represents the best-fit model to the full 
sample shown in Fig. 9. Bottom panel: Same as above, but for the cut 
MegaZ-LRG sample split into the two photometric redshift bins, where 
results for z < 0.529 are shown in black, and for z > 0.529 in red. 
Dotted lines again indicate the best-fit model for the full MegaZ-LRG 
samples, respectively. Note that the red points have been slightly offset 
horizontally for clarity, and that the error bars are correlated. Only the 
data points outside the grey region have been used for the fits. 



exclude more and more of the red sequence. In contrast, for this 
paper, we want to keep all of the red sequence without regard for 
matching the luminosity ranges (indeed, we would like to study 
samples on a wide luminosity baseline in order to measure the 
variation of intrinsic alignments with luminosity). Thus, we wish 
to define a minimum 2 {g - i) that corresponds roughly to that 
for the SDSS spectroscopic LRGs and our revised definition of 
the SDSS Main L3 and L4 red samples. The best choice in this 
context appears to be a cut at 02 (g - i) > 1.7, which should en- 
sure consistency with the other samples within the limits of our 
uncertainty in k + e-corrections. 

To illustrate this cut, we present Fig. 12, which shows two- 
dimensional projections of the relationship between redshift, 
colour, and absolute rest-frame magnitude of the samples. As 
shown, they span a wide range of redshifts (0.05 < z < 0.7) and 
of luminosities (four magnitudes), and with the imposition of 
this new colour cut, the colour ranges are quite similar. MegaZ- 
LRG shows the largest scatter to redder colours; however, this is 
expected given that, as the highest redshift sample, they have the 



largest photometric errors which significantly widens the colour 
distribution at the red end where the g band flux is often only 
weakly detected, especially once the 4000A break moves from 
g to r band. The result of this cut is to reduce the MegaZ-LRG 
sample to 70 % of its original size. The typical sample redshift 
does not change, and the mean luminosity increases marginally 
by 2 % compared to the values of the full MegaZ-LRG sample, 
see Table 3. 

We repeat the intrinsic alignment amplitude fits of Sect. 5.3 
for the cut MegaZ-LRG samples, including the contributions by 
galaxy-galaxy lensing and magnification-shear correlations. We 
continue to use the relation between photometric and spectro- 
scopic redshifts from the 2SLAQ verification sample of Sect. 2.4 
because we do not observe any significant effect by the colour 
cut on this relation, e.g. neither the mean nor the scatter of 
the distribution of differences between photometric and spectro- 
scopic redshifts change beyond the lcr level. The resulting cor- 
relation functions with the best-fit models are plotted in Fig. 13, 
and the corresponding best- fit values for A listed in Table 3. 
For comparison we also show the best-fit models to the uncut 
MegaZ-LRG samples in the figure. Since the cut and uncut sam- 
ples have the same redshifts and luminosities to high accuracy, 
we can ascribe any difference in the signals to a dependence of 
intrinsic alignments on galaxy colour. 

For both the high- and low-redshift MegaZ-LRG samples as 
well as the full sample we find an increase in A for the cut sample 
which has higher 2 (g-i), thus being in line with the expectation 
that redder galaxies have stronger intrinsic alignments. The in- 
crease amounts to 1 1 % for the low-redshift sample and 28 % for 
the high-redshift sample (the corresponding error bars feature 
a similar increase), suggesting a stronger colour dependence at 
higher redshift. However, it should be noted that all observed 
changes due to the colour cut in g - i remain within the 1 <x errors 
and are therefore not statistically significant. 

5.5. Intrinsic alignment model fits to combined samples 

Having addressed the question of the compatibility of the sam- 
ples, we repeat the fits to w„+ for different combinations of 
galaxy samples, now allowing for an additional redshift and lu- 
minosity dependence according to the extended model 



P™ de \k,z,L) 



Ab s P sl (k,z) 



1 +z 
1 +Zo 



pother 



(19) 



where zo — 0.3 is an arbitrary pivot redshift, and Lq is the pivot 
luminosity corresponding to an absolute r band magnitude of 
-22 (passively evolved to z — 0). The matter-intrinsic power 
spectrum Psi is given by the NLA model with the modified red- 
shift dependence discussed in Appendix B, including the nor- 
malisation to SuperCOSMOS. This model contains three free 
parameters {A,/3, Mother), an d, as before, has a fixed dependence 
on transverse scales. 

The amplitude parameter A and the luminosity term can be 
taken out of all integrations leading to w g+ because they neither 
depend on redshift nor comoving distance, so that the parameters 
A and /3 can be varied in the likelihood analysis with low com- 
putational cost. The extra redshift term containing 77 ot her depends 
on the integration variable in (5) though. To facilitate the like- 
lihood analysis for the MegaZ-LRG samples with photometric 
redshifts, we assume that this term can be taken out of the inte- 
gration and is evaluated at the mean z m = (ft + z.2)/2 of the two 
redshifts entering (5). This approximation holds to fair accuracy 
because the corresponding redshift probability distributions are 
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Fig. 14. Constraints by the joint fit to MegaZ-LRG and SDSS spectroscopic samples on the amplitude A of the intrinsic alignment model, the 
extra redshift dependence with power-law index pother, and the index /3 of the luminosity dependence. In the lower left panels, the two-dimensional 
lcr and 2tr confidence contours are given, marginalised in each case over the parameter not shown with flat priors in the range A e [0;20], 
Mother e [-10 : 10], and /? 6 [-5; 5]. The upper right panels display the constraints on A, Mother, an d eac h marginalised over the two remaining 
parameters. Red lines are obtained via fits to the six SDSS LRG samples, green lines for the SDSS LRG and Main samples, yellow lines via fits 
to the MegaZ-LRG and SDSS Main samples, blue lines for MegaZ-LRG and SDSS LRGs combined, and black lines result for the joint fit to the 
MegaZ-LRG, SDSS LRG, L4, and L3 samples. 



sufficiently narrow, given the small photometric redshift uncer- 
tainty for the MegaZ-LRG sample. 

The additional redshift dependence is then integrated over 
in the averaging process in (12) and (14) for photometric and 
spectroscopic samples, respectively. As also the luminosity dis- 
tributions of the galaxy samples under consideration are com- 
pact and narrow, we use the mean luminosity in the luminosity 
term in (19) instead of integrating {L/L^f over the full distribu- 
tion. This is a good approximation even for the full MegaZ-LRG 
sample, which features the broadest luminosity distribution, the 
deviation being below 2 % close to the best-fit values for /3 that 
we determine below. 

We consider joint fits to several combinations of the six 
SDSS LRG subsamples, the two MegaZ-LRG low- and high- 
redshift samples with the colour cut, and the two SDSS Main 



L4 and L3 samples. The resulting two-dimensional marginalised 
confidence contours and marginal one-dimensional posterior 
distributions for the parameter set {A,/3, pother) are shown in 
Fig. 14. The corresponding marginal lcr errors on these parame- 
ters and the goodness of fit are given in Table 4. In the computa- 
tion of marginalised constraints we assumed by default flat priors 
in the ranges A e [0; 20], ^ther e [-10 : 10], and f3 e [-5 : 5]. 
For the combination of the MegaZ-LRG and SDSS Main sam- 
ples, which yields weak and degenerate constraints, we extend 
the prior range of pother to +20. Note that in this case the poste- 
rior has not yet decreased very close to zero at/3 = 5, but still we 
expect the influence of the /3-prior on the marginal constraints to 
be negligible. 

Combining all SDSS LRG samples we can constrain /3 well, 
i.e. the power-law slope of the luminosity evolution of the in- 
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Table 4. lcr marginalised confidence limits on the amplitude A, the slope of the additional redshift dependence ^ ot h e r, and the luminosity dependence 
/3 of the intrinsic alignment model (19). 



sample combination 


*V sample 


A 


Mother 


P 




p(>x z ) 


SDSS LRG 
SDSS LRG + L4 + L3 
MegaZ-LRG + L4 + L3 
MegaZ-LRG + SDSS LRG 
MegaZ-LRG + SDSS LRG + L4 + L3 


6 

8 
4 
8 

10 


6 27 +110 
6 32+ 105 

^ 09+0:77 
J y -0.75 
5 76+ - 60 


9 97+2.07 
^■^'-2.06 
9 97+I.93 

-6-80 + ?^ 

"•^'-0.96 

-0 27 +aso 

"""-0.79 


1 17+O.30 
1-1 -0.32 

2 87 + ^' 
1 10+^ 

1 - lu -0.30 

1 13 +0 - 25 

1-1 -0.20 


1.37 
1.23 
0.60 
1.15 
1.09 


0.09 
0.16 
0.89 
0.24 
0.32 



Notes. Shown are the fit results for A, pother, and f} for various sample combinations. In addition we list the number of subsamples iV sam pie use d 
for the joint fit and the reduced x 2 including p-value. The number of the degrees of freedom is given by 5 7V siimple - 3. We have used flat priors 
A e [0; 20], t] oAs . t e [-10; 10], and /? e [-5; 5], except for the case 'MegaZ-LRG + L4 + L3\ where 77 oCllcr 6 [-20; 20] was assumed instead. 



trinsic alignment amplitude, while the errors on 77 0t her remain 
large, as expected for the short baseline in redshift. Adding in 
the two redshift-binned MegaZ-LRG samples greatly improves 
constraints on the extra redshift evolution and also narrows down 
the possible values of A. The marginalised lcr contours for the 
SDSS LRG only and the SDSS + MegaZ-LRG fits are consis- 
tent, although the MegaZ-LRG data clearly prefers smaller val- 
ues of pother- The reduced^ 2 improves by 16 % when adding the 
MegaZ-LRG data. Incorporating in addition the SDSS Main L3 
and L4 samples at low redshift and with low mean luminosity 
further tightens constraints on all parameters, in particular the 
amplitude A. The SDSS Main samples have small constraining 
power due to the large error bars on their signals. Their inclu- 
sion yields largely consistent confidence regions in combination 
with both SDSS LRG and MegaZ-LRG samples and decreases 
the reduced^- 2 by another 5 %. 

It is not a priori clear that the intrinsic alignment model deter- 
mined for the LRG samples also holds for the fainter, non-LRG 
SDSS Main samples. The validity of this model for galaxies with 
luminosities around and below L* is paramount for its applica- 
bility to intrinsic alignments in cosmic shear data that has many 
more faint red galaxies than LRGs. To affirm consistency, we 
read off the combination of intrinsic alignment parameters that 
yield the minimum x 2 f° r the joint fit to the SDSS LRG and 
MegaZ-LRG samples. Then we compare the^- 2 of this parame- 
ter combination for the fit to the combined L3 and L4 samples to 
the minimum^ 2 obtained for fitting to this sample combination. 
We find a difference A^- 2 = 2.01 to the minimum^- 2 of the latter 
samples in the parameter range A e [0; 20], 77 0t her £ [-10 : 10], 
and /3 e [-5; 5], corresponding to a p-value of 0.57. Thus the 
SDSS Main samples are fully consistent with the LRG data sets, 
but note that the L4 and L3 data are noisy and hence not very 
conclusive, e.g. we also compute a reduced^- 2 of 1.19 (L3) and 
0.89 (L4) for a fit to zero. For illustration, we have plotted both 
the best-fit models for a fit to L3 and L4 only, and for the fit to 
the SDSS LRG and MegaZ-LRG samples in Fig. 15. 

The joint fit of all considered samples clearly favours an in- 
crease of the intrinsic alignment signal with galaxy luminosity; 
indeed we find that /3 < 0.5 can be excluded at more than the 
Act level. The data are perfectly consistent with 77 ot h er = 0, i.e. 
with the redshift evolution inherent to the corrected NLA model 
as discussed in Appendix B. The combination of MegaZ-LRG 
and SDSS LRG samples is the main driver for these constraints 
on the redshift dependence. Using the NLA model with the red- 
shift dependence of Hirata & Seljak (2004) instead would have 
shifted the constraints on 77 ther by -2. Consequently, our results 
disfavour this redshift dependence of the original NLA model, 
excluding pother > 2 by about 2.8cr. We find an overall intrin- 




r p [Mpc/h] 

Fig. 15. Consistency between the combined LRG and SDSS Main sam- 
ples. Shown is w g+ for the red SDSS L4 (in black) and L3 (in red) sam- 
ples. The solid curves correspond to the best-fit models to the com- 
bined L3 and L4 samples when varying A, 77 ot her, and /?. The corre- 
sponding best-fit models to the combined SDSS LRG and MegaZ-LRG 
samples are given by the dotted lines. The prior ranges A e [0;20], 
Mother s [-10; 10], and /3 6 [-5; 5] were employed in both cases. Note 
that the red points have been slightly offset horizontally for clarity, and 
that the error bars are correlated. 

sic alignment normalisation A = 5.8 + 0.6, which e.g. trans- 
lates into an amplitude of 90 % of the standard NLA model with 
corrected redshift dependence and SuperCOSMOS normalisa- 
tion for a typical red galaxy with L = 0.2 Lq « L, at redshift 
z = 0.5 (see Appendix C for a justification of these values). 
Using the NLA model with the redshift dependence of Hirata 
& Seljak (2004), which has been employed in most weak lens- 
ing studies hitherto, yields an amplitude of about 40 % of the 
SuperCOSMOS normalisation. 

Our findings for the marginalised parameter constraints on 
/3 from the fits to the SDSS LRG samples are in good agree- 
ment with the results presented in Hirata et al. (2007), despite 
a different binning in transverse separation and a fit of a pure 
power-law dependence on z and r p instead. The latter yields a 
scaling with r p which is comparable to the NLA model. Due to 
these differences in modelling, the intrinsic alignment amplitude 
parameters are not directly compatible, while the Hirata et al. 
(2007) power-law index for the redshift term should roughly cor- 
respond to 77 0t her _ 1, see (B.7). The results for the redshift de- 
pendence are also consistent at the lcr level, albeit with large 
error bars. In addition Hirata et al. (2007) considered joint con- 
straints from SDSS LRG and 2SLAQ samples, but did not im- 
pose colour cuts on 2SLAQ, so that a quantitative comparison 
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Fig. 16. k + e-corrections to z = using different sets of templates. 
Top panel: k + e-corrections in the colour g - i used to determine the 
compatibility in colour ranges of the galaxy samples. The black curve 
corresponds to the templates used in this work by Wake et al. (2006) 
which are based on the stellar population synthesis code by Bruzual & 
Chariot (2003) using passive evolution. If a low level of continuous star 
formation is added to the model, the red solid line results. In addition we 
show the k + e-corrections from Banerji et al. (2010) using population 
synthesis models by Maraston et al. (2009) as the blue dotted curve. 
Bottom panel: k + e-corrections in the r band in which we compute 
galaxy luminosities. The coding of the curves is the same as above. 
The grey vertical lines indicate typical redshifts of the samples under 
consideration. 



is difficult. However, we observe similar tendencies in the best- 
fit values of the intrinsic alignment parameters when adding the 
2SLAQ and MegaZ-LRG samples, respectively. Besides, param- 
eter errors from the joint analysis of SDSS LRG and 2SLAQ, or 
SDSS LRG and MegaZ-LRG samples are of the same order of 
magnitude. 

5.6. Dependence on k + e -corrections 

We have specified our intrinsic alignment model in terms of the 
physically meaningful rest-frame r band luminosity (passively 
evolved to z — 0), which renders our fit results dependent on 
the k + e-corrections employed. Recently, Banerji et al. (2010) 
have utilised new corrections based on improved stellar popula- 
tion synthesis models (Maraston et al. 2009) that perform better 
than those of Wake et al. (2006) at describing observed colours 
of LRGs. We compare these different k + e-corrections in the 
r band in Fig. 16, finding significant differences in particular at 
high redshift. 



The difference between the results from Wake et al. (2006) 
and Banerji et al. (2010) has a negligible effect on the luminosi- 
ties computed for the SDSS Main samples at low redshift, but 
amounts to approximately 0.1 mag for the SDSS LRG samples 
and, the details depending on the population synthesis model 
used, to about 0.2 mag at the mean redshift of the MegaZ-LRG 
sample. Hence, switching to luminosities calculated according to 
Banerji et al. (2010) would imply a 10 % increase in the intrinsic 
alignment amplitude A for the fits to the combined SDSS LRG 
samples, using the corresponding best-fit value for f3. While this 
change is subdominant to the Icr error, a change by 0.2 mag for 
the MegaZ-LRG sample leads to a shift by about 20 % in the 
term (L/Lof, which clearly exceeds the lcr errors e.g. of the 
joint fit by all galaxy samples. Besides, the effect is redshift- 
dependent, so that all three fit parameters A, fl, and Mother would 
be affected in that case. 

However, in principle the definition of luminosity used for 
(19) is arbitrary, as long as it is consistent for all samples. This is 
the case since we employ the A; -He-corrections according to Wake 
et al. (2006) throughout. As demonstrated above, the transforma- 
tion of our intrinsic alignment model to a different convention 
for galaxy luminosity has a considerable effect on all parameters 
and needs to be executed with care. 

The actual quantity of interest is not the power spectrum (19) 
but rather the observable intrinsic alignment signals, in particular 
those contaminating cosmic shear surveys. Changes in the val- 
ues of the intrinsic alignment fit parameters are only meaningful 
in how they modify these observables. The observable intrin- 
sic alignment signals are determined by observer-frame apparent 
magnitude limits of a survey and can consistently be calculated 
from (19) if the same k + e-corrections as used in the computa- 
tion of galaxy luminosities for the model are applied. In Sect. 6 
we will use such a procedure to predict the contamination of a 
cosmic shear survey by intrinsic alignment signals derived from 
our best-fit model. As will be detailed in there, there are sources 
of uncertainty, especially concerning the choice of luminosity 
functions, that are likely to be more important than the effect of 
k + e-corrections. 

Besides the luminosities of the galaxy samples, their 2 (g-i) 
colours used to assess the compatibility of the samples are ex- 
pected to depend on the k + e-corrections as well. Again, we 
have made sure that identical k + e-corrections were employed 
to produce Fig. 12, so that the MegaZ-LRG samples are consis- 
tent with the other data also in their colour cuts. Note however 
that, contrary to the luminosity dependence, our results cannot 
readily be reformulated for galaxy colours obtained from other 
versions of k + e-corrections due to the imposed colour cut which 
is dependent on the Wake et al. (2006) templates. 

As shown in Fig. 16, we have compared the differences be- 
tween the k + e-corrections in the g and i bands obtained from 
Wake et al. (2006) and Maraston et al. (2009, see Banerji et al. 
2010 for details) as a function of redshift. We find that the latter 
are in good agreement with the two variants of colour correc- 
tions derived from Wake et al. (2006). The Maraston et al. (2009) 
k + e-corrections display a slightly steeper increase with redshift, 
yielding g - i corrections that are 0.03 lower at the redshift of the 
SDSS Main samples and about 0. 1 higher at the mean redshift of 
the full MegaZ-LRG sample compared to the mean of the Wake 
et al. (2006) models. Consequently, the MegaZ-LRG samples 
would on average be bluer with respect to the other samples if 
we had used the k + e-corrections based on Maraston et al. (2009) 
instead. While a difference of 0. 1 in the colour appears to be sub- 
stantial compared to the typical spread of the MegaZ-LRG sam- 
ple in 2 (g - i) of > 0.5, the uncertainty in the colour cut due to 
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the fuzziness in the lower g - i limit of the SDSS LRG and Main 
samples is of the same order, see Fig. 12. Besides, as Fig. 16 top 
panel suggests, this level of uncertainty due to the differences 
between the templates by Wake et al. (2006) and Maraston et al. 
(2009) is of the same order as the uncertainty in the evolutionary 
model chosen for a galaxy within a given set of templates. 



5.7. Systematics tests 

As in Mandelbaum et al. (2006), Hirata et al. (2007), and 
Mandelbaum et al. (2010) we perform several tests to ensure 
that our results from the MegaZ-LRG data are not contaminated 
by instrumental or other effects. We also repeat the systematics 
tests for the re-defined red SDSS Main L3 and L4 samples. No 
signatures of systematics were found during the previous analy- 
ses of the SDSS LRG and the original SDSS Main samples (see 
the aforementioned references). 

First, we compute the correlation function w gx using the 
cross-component of the shear instead of the radial component. 
The cross-component and thus w gx change sign under parity 
transformations, measuring a net curl of the galaxy shape dis- 
tribution. Since we do not expect that galaxy formation and evo- 



e 0.5 



X 



lution violates parity symmetry, any non-vanishing w gx serves * 
as an indicator for systematics, such as residual PSF distortions. j 
The resulting signals for the full MegaZ-LRG sample, as well 
as for the two MegaZ-LRG redshift bins, are shown in the cen- 
tre panel of Fig. 17. All signals are comfortably consistent with 
zero, and fits of a zero line to the data in the full range of r p 
yield reduced x 1 values well below unity and correspondingly 
large p-values, as shown in Table 5. We repeat this analysis for 
the same MegaZ-LRG samples but with the cut in g - i imposed, 
again finding no evidence for systematics, see the bottom panel 
of Fig. 17. 

Furthermore we consider w g+ computed only for large line- 
of-sight separations n at which one does not expect astrophys- 
ical correlations anymore. A non-zero signal in this measure 
could for instance be caused by an artificial alignment of galaxy 
images in the telescope optics. Due to the photometric redshift 
scatter in MegaZ-LRG data, we use much larger values of | XT | 
for this systematics test than preceding works, integrating the 
three-dimensional correlation function along the line of sight 
for 270 < ini/t/r'Mpc] < 315. Still, gl correlations might not 
be completely negligible. We estimate the signal from the best- 
fit intrinsic alignment model obtained in the foregoing section, 
finding amplitudes below 0.004 for r p > 6/z~ 1 Mpc and below 
0.018 for all r p scales considered. Thus, residual gl correlations 
should be negligible in this case, and indeed the resulting w g+ is 
consistent with zero. Note that we cannot apply this test to w gg 
because the much stronger galaxy clustering signal is not negli- 
gible even in the extreme range of n that we have chosen. 

We also compute w gx , as well as w g+ in the range 100 < 
\Tl\/[h- l Mpc] < 150, for the SDSS L3 and L4 samples with the 
new colour cut to isolate red galaxies, as shown in the top panel 
of Fig. 17. Note that for these spectroscopic samples we can re- 
sort to much smaller values of |FI| to obtain a range which is not 
expected to contain physical correlations anymore. All system- 
atics tests for the SDSS Main L3 and L4 samples are consistent 
with zero, see Table 5. Although some points deviate from zero 
clearly outside the lcr limit, and reduced^ 2 values slightly ex- 
ceed unity, the p-values indicate that there is no significant signal 
in the data. 
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Fig. 17. Projected correlation functions w gx and w g+ for large line-of- 
sight separations as a function of comoving transverse separation r p . 
Top panel: Shown is w gx for the re-defined SDSS Main L4 (red points) 
and L3 (black points) samples. Moreover w g+ integrated along the line 
of sight for 100 < |n|/[/i _I Mpc] < 150 is plotted as grey points for 
L3 and orange points for L4. Note that the orange, black, and red 
points have been slightly offset horizontally for clarity. Centre panel: 
Shown is w gx for the full MegaZ-LRG sample in grey, for the low- 
redshift sample with z < 0.529 in black, and for the high-redshift sam- 
ple with z > 0.529 in red. In addition we plot the correlation function 
Wg + for the full MegaZ-LRG sample, integrated along the line of sight 
for 270 < |n|/[/2~'Mpc] < 315, as green points. Note that the black, 
red, and green points have been slightly offset horizontally for clarity. 
Bottom panel: Same as in centre panel, but for the MegaZ-LRG sample 
with the colour cut imposed. Note that we did not repeat the computa- 
tion of w g+ for large \Yi\. Error bars at different r p are correlated. Note 
that all correlation functions have been rescaled by [r,,/(l /i~'Mpc)] 06 
for easier inspection. All signals are consistent with zero. 



As expected, we find that error bars on w gx are of similar 
size as for w„+ when calculated for the same range in line-of- 



sight separation; compare e.g. the centre panel of Fig. 17 to Fig. 9 
(note the different scaling of the ordinate axes). Since the pro- 
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Table 5. Systematics tests for the MegaZ-LRG data and the re-defined SDSS L3 and L4 samples. 



sample 


correlation 




' v red v / 




y 1 ( 


r p > 6Mpc/A) 


p(>x z ) 


MegaZ-LKLr 


Wg X , all z, (yu) 




n ^ 1 
U.3 1 


U.oS 




a no 




MegaZ-LKLr 


... ^ ^ n con 

Wgx, z < 0.529, 


(90) 


i no 
1.08 


0.3 / 




0.50 


0. / 1 


Mera7-T RG 


w 7 > 5?9 




0.38 


0.96 




0.49 


0.79 


MegaZ-LRG 


w g +, all z, large 


mi 


1.01 


0.43 




1.44 


0.20 


MegaZ-LRG (g - i cut) 


Wgx, all z, (90) 




0.72 


0.70 




0.77 


0.57 


MegaZ-LRG (g - i cut) 


Wgx, z < 0.529, 


(90) 


1.25 


0.25 




1.66 


0.14 


MegaZ-LRG (g - i cut) 


w fc ,x, z > 0.529, 


(90) 


0.30 


0.98 




0.22 


0.96 


SDSS Main L3 red 


Wgx 




0.97 


0.47 




0.70 


0.62 


SDSS Main L3 red 


w g+ , large |II| 




0.49 


0.90 




0.21 


0.96 


SDSS Main L4 red 


Wgx 




1.24 


0.26 




0.49 


0.78 


SDSS Main L4 red 


w g+ , large |II| 




1.06 


0.39 




1.01 


0.41 



Notes. We list the reduced^ 2 and the corresponding p-values of a zero signal fit to the correlation function w gx , and to w g+ for large values of 
|II|. We use 270 < ini/r/r'Mpc] < 315 for the photometric MegaZ-LRG data and 100 < m/ITr'Mpc] < 150 for the spectroscopic SDSS Main 
samples. We include all transverse scales available into the fit (third and fourth columns), but also restrict r p to scales larger than 6/r'Mpc (fifth 
and sixth columns), i.e. the range in which the intrinsic alignment fits were performed. The numbers in parentheses indicate IT max in units of 
/r'Mpc for the MegaZ-LRG samples. 



jected correlation functions are given by the sum of bins in II, 
and not the mean, see e.g. (12), one expects larger errors on w gx 
or w g+ if the range in II included is broader. Therefore w g+ in 
the range 100 < m/t/r'Mpc] < 150 for the SDSS Main sam- 
ples and w g+ in the range 270 < |Il|/[/z _1 Mpc] < 315 for the 
MegaZ-LRG sample features smaller error bars than w gx , calcu- 
lated in the range ini/f/T^pc] < 60 and m/t^Mpc] < 90, 
respectively. 

One additional type of systematic is the calibration of the 
shear. A multiplicative calibration offset would manifest directly 
as a multiplicative factor in w g+ . Aside from any impact on the 
best-fitting intrinsic alignment amplitude A from the combined 
samples, there could also be some effect that is a function of 
galaxy properties (typically apparent size and S /N, see Massey 
et al. 2007a), which would manifest as a difference between the 
SDSS Main samples (very bright apparent magnitudes and well- 
resolved), the spectroscopic LRGs (moderately bright and well- 
resolved), and MegaZ-LRG (significantly fainter and less well- 
resolved). Because the samples occupy different places in both 
redshift and luminosity space, the effect of such a bias cannot 
be estimated in a straightforward way. However, the shape mea- 
surements used for this work were subjected to significant sys- 
tematics tests in Mandelbaum et al. (2005), including tests for 
calibration offsets between different samples, and thus we do not 
anticipate that shear calibration is a significant systematic rela- 
tive to others (such as photometric redshift error uncertainties, 
or k + e-corrections) or relative to the size of the statistical error 
bars on w e+ measurements. 



6. Constraints on intrinsic alignment contamination 
of cosmic shear surveys 

Intrinsic alignments constitute the potentially major astrophysi- 
cal source of systematic uncertainties for cosmic shear surveys. 
If left untreated, they can severely bias cosmological parameters 
estimates (e.g. Bridle & King 2007). If the contamination by in- 
trinsic alignments is well known, it can ideally be incorporated 
into the modelling by subtracting the mean intrinsic alignment 
signal from the lensing term and accounting for the residual un- 
certainty in the systematic by introducing nuisance parameters 
over which one can then marginalise. 



To elucidate the implications for cosmological constraints 
from cosmic shear surveys by our constraints on intrinsic align- 
ments, we will optimistically assume that the mean systematic 
signal is indeed given by our best-fit model. In this approach, the 
decisive quantity is not the mean value of the bias on cosmolog- 
ical parameters, which can be easily corrected for by subtracting 
the mean intrinsic alignment signal, but the uncertainty in the 
bias due to uncertainty in intrinsic alignment model parameters, 
which directly affects the accuracy with which the cosmologi- 
cal model can be constrained when taking the systematics into 
account. We do not address uncertainty due to adoption of the 
generalised NLA model (19), i.e. the possibility that the under- 
lying intrinsic alignments model is different, because we see no 
tension between the model and our data. 

We assess the range of possible biases on cosmological pa- 
rameters that originate from intrinsic alignment signals using the 
constraints obtained from the foregoing investigation. Since the 
SDSS Main L3 and L4 samples proved to be fully consistent 
with the results for the two LRG samples, we will assume in the 
following that our intrinsic alignment model also holds for typi- 
cal, less luminous early-type galaxies predominantly found in a 
cosmic shear survey. We emphasise that this study requires the 
extrapolation of the best-fit intrinsic alignment model to combi- 
nations of galaxy redshifts and luminosities that have not been 
probed directly by any of the galaxy samples analysed in this 
work. However, that extrapolation is less worrisome now that 
we have galaxy samples at z ~ 0.6, given that the lower redshift 
galaxies in a cosmic shear survey will tend to be the greatest cul- 
prits in causing GI contamination due to the scaling of the effect 
with redshift separations of galaxy pairs. 

By means of a Fisher matrix analysis we compute the ef- 
fect on a present-day, fully tomographic (i.e. including all in- 
dependent combinations of redshift bins) cosmic shear survey, 
roughly following CFHTLS parameters (Hoekstra et al. 2006; 
Semboloni et al. 2006; Fu et al. 2008). To calculate the matter 
power spectrum, we use the same cosmology, transfer function, 
and nonlinear correction as outlined in Sect. 4. 1 . For computa- 
tional simplicity we use the convergence power spectrum (the 
GG signal henceforth) as the observable cosmic shear two-point 
statistic, using the following Limber equation, 

^.^rf^M,,^). (20) 



23 



B. Joachimi et al.: Constraints on intrinsic alignment contamination of weak lensing surveys using the MegaZ-LRG sample 



where again denotes the lensing weight. Instead of specify- 
ing a photometric redshift for a redshift probability distribution, 
we switch here to the usual notation of using an index i that 
characterises a (broad) distribution /? w (z) entering (9). The cor- 
responding Limber equations for the GI and II signals can be 
readily formulated accordingly, yielding (e.g. Hirata & Seljak 
2004) 
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where we assume that the intrinsic shear power spectrum can 
be described by an extension of our intrinsic alignment model 
analogous to (19), 

l2^„,her / j \2P 
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(k,z,L)=A z P u (k,z) 



1 + z 

1 +ZQ 



(23) 



where P u (k,z) is given by the NLA model, i.e. (B.6) with the 
linear matter power spectrum replaced by the full, nonlinear one. 
Writing (23) with the square of the extra redshift and luminosity 
terms (again similar to Kirk et al. 2010) involves the assumption 
that the galaxies correlated are located at similar redshifts, which 
is valid because the II signal is restricted to physically close 
pairs; see also the narrow kernel in (21). Then the two galax- 
ies of a pair also underlie the same luminosity distribution, and 
since we average (23) over this distribution, one can write lP 
to good approximation although the luminosities of the galax- 
ies in an individual pair may be largely different. Note that we 
have not measured intrinsic ellipticity correlations in this work, 
so that the computation of the II signal is entirely based on the 
validity of the NLA model. However, the contribution of the II 
signal to the bias on cosmology will be smaller than the one by 
GI correlations in a CFHTLS-like survey, so that this assumption 
does not affect our conclusions substantially. 

We employ an overall redshift distribution according to 
Smailetal. (1994), 



Ptot(z) <* 
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(24) 



with parameters a Smai i = 0.836, /? Sma ii = 3.425, and z S maii = 
1.171 yielding a median redshift of 0.78 (Benjamin et al. 2007). 
We slice this distribution into 10 'photometric' redshift bins such 
that every bin contains the same number of galaxies. The corre- 
sponding redshift distribution for each bin is then computed via 
the formalism detailed in Joachimi & Schneider (2009), assum- 
ing a Gaussian photometric redshift scatter of width 0.05(1 + z) 
around every spectroscopic redshift. We compute Gaussian co- 
variances for the power spectra including cosmic variance and 
shape noise (for details see Joachimi et al. 2008), assuming a 
survey size of A survey = lOOdeg 2 . Shape noise is incorporated 
with an overall galaxy number density of «n = 12 arcmin~ 2 and 
a dispersion of the absolute value of the complex intrinsic ellip- 
ticity of 0.3 (Hoekstra et al. 2006). 

We consider a parameter vector p = {Q. m , cr g , h, n s , Qb, wq) 
for the cosmological analysis, for a flat universe with constant 
dark energy equation-of-state parameter wq. Assuming that the 
covariance is not dependent on these parameters (see Eifler et al. 
2009), one obtains the Fisher matrix (Tegmark et al. 1997) 

^ = ^ ^ C ov-(C^),0.))^, (25) 



where we use 40 logarithmically spaced angular frequency bins 
between ( — 10 and ( = 3000. With the Fisher matrix, one 
can calculate the bias on a cosmological parameter via (e.g. 
Kim et al. 2004; Huterer & Takada 2005; Huterer et al. 2006; 
Taylor et al. 2007; Amara & Refregier 2008; Kitching et al. 
2009; Joachimi & Schneider 2009) 



(26) 



C,i< j,k<l 



x Cov- 1 (d®{€), OO) 



dp v 



l,i< j,k<l 



dp v 



where the systematic is given by the sum of II and GI power 
spectra. Note that the parameter bias is independent of the sur- 
vey size while the statistical errors obtained from F MV are pro- 
portional to l/-y/A survey . 

The intrinsic alignment analysis presented above only dealt 
with red galaxies, whereas a typical galaxy population in cos- 
mic shear surveys is dominated by blue galaxies for which 
Mandelbaum et al. (2010) reported a null detection for a galaxy 
sample spanning a similar range of redshifts as in this paper. 
Thus we assume that only the red fraction f r of galaxies in the 
survey carries an intrinsic alignment signal. Consequently the II 
power spectrum is multiplied by a factor / r 2 , and the GI power 
spectrum by f r , resulting in the same model as used by Kirk et al. 
(2010). Note that this approach is overly simplistic in splitting 
the galaxy population into two disjoint groups with largely dif- 
ferent intrinsic alignment properties although one expects the in- 
trinsic alignment parameters to vary in a more continuous man- 
ner with galaxy colour. 

Moreover it is important to note that we ignore any uncer- 
tainty in the null measurement of blue galaxy intrinsic alignment 
which would add to the scatter of systematic errors on the cos- 
mological parameters. In principle it should be feasible to take 
into account finite constraints on the intrinsic alignment param- 
eters determined from blue galaxy samples, such as those stud- 
ied by Mandelbaum et al. (2010). However, before incorporating 
them into this formalism, these samples would have to undergo 
the same compatibility tests as performed in this work for red 
galaxy samples, and then be combined to yield joint fits. These 
steps are beyond the scope of this analysis and left to future in- 
vestigation. 

We expect that both f r and the distribution of luminosities 
of red galaxies depend on redshift and thus have different val- 
ues in each photometric redshift bin of our fictitious cosmic 
shear survey. To estimate realistic values for these parameters, 
we make use of the luminosity functions provided by Faber et al. 
(2007). They fit Schechter functions <f>(L,z) jointly to samples 
from SDSS, 2dF, COMBO- 17, and DEEP2 in redshift bins out 
to z ~ 1, considering early- and late-type galaxies individually. 
The criteria used by Faber et al. (2007) to separate red and blue 
galaxies differ from the ones employed in this work, but still we 
consider their samples as representative for differentiating be- 
tween blue galaxies with negligible intrinsic alignments and red 
galaxies with an intrinsic alignment signal consistent with our 
best-fit model. We defer the technical details and also further 
discussion of this approach to Appendix C. 

In combination with a minimum galaxy luminosity 
Lmin(z, Him), computed from the apparent magnitude limit at each 
redshift, a set of luminosity functions also specifies the redshift 
probability distribution p to t(z)- However, the luminosity func- 
tions by Faber et al. (2007) are unlikely to reproduce exactly the 
redshift probability distribution of CFHTLS because weak lens- 
ing surveys have additional galaxy angular size cuts and thus 
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deviate from a purely flux-limited sample. Also, we must ex- 
trapolate to fainter magnitudes than the sample used to deter- 
mine the luminosity function, which can be a source of signifi- 
cant uncertainty, particularly for blue galaxies given their steep 
faint-end slope. Besides, the blue galaxy sample of Faber et al. 
(2007) is composed of several galaxy types, so that it is un- 
clear which fe-corrections and filter conversions are applicable. 
Hence, we make the assumption that the red galaxy luminos- 
ity functions from Faber et al. (2007) are compatible with the 
selection criteria of a weak lensing survey and derive the total 
comoving volume density of galaxies from (24) by noting that 
«v = dA7dV com = dA7dz(dV com /dz) _1 . Then the fraction of red 
galaxies reads 

»V,redfc Him) X 2 (z)x'(z) f°° 

Mz,n im ) = — — = — AL c/>(L,z) , (27) 

«v,tot(z) nnpiotiz) J Lm uz,n im ) 

where the lower limit of the integration L m [ n is determined by the 
magnitude limit of the survey which we assume to be r\{ m = 25, 
roughly compatible with CFHTLS (;'ii m = 24.5; Hoekstra et al. 
2006). The redshift dependence of L m [ n is introduced by the con- 
version of r\ lm to absolute magnitude and by the k + e-correction. 
The latter is again computed in the r band by means of the tem- 
plates used in Wake et al. (2006) for early-type galaxies (specif- 
ically, by a ^-correction that is between that of the two very 
closely related templates used in that paper, see also Fig. 16). 

In contrast to the various SDSS samples in which we mea- 
sured intrinsic alignments, the distribution of galaxy luminosi- 
ties in each photometric redshift bin of our mock cosmic shear 
survey is wide, so that our approximation that we can replace L 
by the mean luminosity in (19) and (23) is not sufficiently accu- 
rate anymore, in particular in those regions of the intrinsic align- 
ment parameter space where ft deviates significantly from unity. 
Instead, we have to average (19) and (23) over the luminosity 
function, which reduces to evaluating 

where we have taken into account that the rest-frame reference 
luminosity Lq must be e-corrected back to redshift z, which we 
do via the redshift dependence of M* given in Faber et al. (2007); 
see again Appendix C for details. For every photometric redshift 
bin of the cosmic shear survey, we use the values of ((L/Lo'f'j 
and f r at the median redshift of the underlying redshift distri- 
butions, which is a good assumption if the redshift distributions 
corresponding to the photometric redshift bins are sufficiently 
narrow. 

We further assess the accuracy and limitations of this ansatz 
in Appendix C by comparing results from sets of luminosity 
functions other than those from Faber et al. (2007). We also pro- 
vide volume densities n v , re d and mean luminosities for a range 
of redshifts and limiting r band magnitudes, which can be em- 
ployed together with our best-fit intrinsic alignment model pa- 
rameters to estimate the expected intrinsic alignment contami- 
nation of other surveys. 

The GI and II signals are then computed via (19) and (23), 
where the free intrinsic alignment parameters {A,/3, pother) are de- 
termined as follows. We overlay the three-dimensional lcr and 
2cr volumes of the intrinsic alignment fits to four of the sample 
combinations shown in Table 4 and Fig. 14 with a square grid, 
containing N nodes in total. For the combination of \A,/3, pother) 
on each grid node we compute the projected intrinsic alignment 
power spectra according to (21) and (22) and subsequently the 



parameter biases via (26). This way we obtain a bias vector 
b = \b(p\), .. ,b(pjf D )} T , where is the number of cosmolog- 
ical parameters under consideration, in cosmological parameter 
space for every grid node within the 1 cr or 2<x confidence volume 
in intrinsic alignment parameter space. 

We convert the ensemble of N parameter bias vectors 
{Ai, .. into a distribution of bias values via Gaussian ker- 
nel density estimation, i.e. we approximate this distribution by 

1 ^ 1 f (xj-biQfj)) 2 } 

where we use No = 2 when considering the distribution in a two- 
dimensional parameter plane, and = 1 when computing the 
one-dimensional distributions. The widths A of the Gaussians 
in every dimension of cosmological parameter space are free pa- 
rameters, and we choose them to take the minimum values which 
still produce a smooth distribution (except for small wiggles in 
some of the sparsely sampled regions). While we use six cos- 
mological parameters to compute the biases on cosmology, we 
focus in our presentation of the uncertainty in the biases on a 
subset with three cosmological parameters of particular inter- 
est in cosmic shear analyses, {£2 m , cr 8 , wq}. For the tightly con- 
strained parameters £2 m and cr$ we use A = 0.001 and in the 
dimension corresponding to wq we set A = 0.005. Note that we 
use the same widths for all sample combinations considered in 
order not to distort the comparison between the resulting bias 
distributions. 

In Fig. 18 we show the contours comprising 99 % of the dis- 
tribution (29) in the two-dimensional parameter planes spanned 
by all pair combinations in the set \Q. m , erg, wq}, sampling from 
the posteriors of the intrinsic alignment parameters obtained for 
SDSS LRGs alone, SDSS LRGs and SDSS Main samples com- 
bined, SDSS LRGs and MegaZ-LRG combined, and the joint 
analysis of MegaZ-LRG, SDSS LRG and SDSS Main samples. 
In this figure we have given the parameter biases (and not the 
parameter values) on the axes such that in the absence of any in- 
trinsic alignment contamination, the contours should be centred 
around (0; 0) in each panel. 

The general direction of parameter biases, for instance along 
the £2 m - cr 8 degeneracy, is in agreement with other predictions 
on biases due to intrinsic alignments (see for instance Joachimi 
& Schneider 2009) 11 . As is evident from (26), if the GI term 
dominates, which is expected for deep cosmic shear surveys, the 
bias is proportional to the amplitude parameter A of the intrinsic 
alignment model. Thus the remaining uncertainty in A explains 
the strong elongation of the contours, pointing approximately 
radially away from (0; 0). The large errors on intrinsic alignment 
parameters, in particular on pother, in the case of using the SDSS 
LRG samples alone allow for a vast region of possible parameter 
biases. Adding the SDSS L4 and L3 samples slightly narrows the 
contours, but does not reduce their radial extent. 

The contours tighten dramatically when adding in the 
MegaZ-LRG data which allowed us to fix the redshift depen- 
dence to good accuracy. The additional information provided by 
the SDSS Main samples constrains the total amplitude of the in- 
trinsic alignment signal still better, thereby further reducing e.g. 

1 1 Note that in deep cosmic shear surveys the negative GI signal domi- 
nates the contamination by intrinsic alignments, so that one may expect 
a negative bias on the power spectrum normalisation <x 8 . However, the 
bias acts mainly along the £l m - <x 8 degeneracy line (apparently causing 
a positive systematic shift in tr 8 ), with only a small shift perpendicular 
to this line in the downward direction. Consequently, for fixed £l m , the 
bias on <x 8 is indeed negative. 
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Fig. 18. Bias on cosmological parameters due to intrinsic alignments for a CFHTLS-like weak lensing survey. Shown are the regions in which 
99% of the possible biases on the parameters {Q. m ,crg,w } are located when the parameters {A,/3, J7 0t her} m the intrinsic alignment model are 
sampled from the l<r confidence region (thick lines) and the 2<r confidence region (thin lines) of our fits. The contours resulting from the SDSS 
LRG constraints are shown in red, the ones from the SDSS LRG + SDSS Main (L3 and L4) constraints in green, the ones from the MegaZ- 
LRG + SDSS LRG constraints in blue, and the contours from the joint constraints by the MegaZ-LRG, SDSS LRG, and SDSS Main samples 
in black. The grey regions indicate the lcr confidence regions of the constraints on cosmological parameters. For this analysis the 6 parameters 
{Q. m ,crg,h, n s ,Sl b , vv ) were varied. The fraction f r of red galaxies and the distribution of luminosities of red galaxies were determined from the 
luminosity functions provided by Faber et al. (2007). Note that the contours corresponding to the 2<r confidence region for the SDSS LRG-only 
and SDSS LRG + SDSS Main constraints extend far beyond the plot boundaries. 



the extent of the 2cr contours by about a factor of two along the 
degeneracy direction. Note that, since the 2<x regions of the in- 
trinsic alignment model fits to SDSS LRGs only and SDSS LRG 
and MegaZ-LRG samples combined do not completely overlap 
(see Fig. 14), the corresponding bias distributions partly cover 
different regions in cosmological parameter space. 

In Table 6 we list the 1 <x marginalised statistical errors for the 
three cosmological parameters of interest, obtained from (25) via 
CstntiPfj) = ^(F^ 1 )^- Moreover we give the size of the interval 
that contains 99 % of the one-dimensional distribution (29) when 
sampling from the lcr and 2<x confidence volumes of the intrin- 
sic alignment parameter fits, again as a measure for the spread 
of biases on cosmology. In agreement with the two-dimensional 
plots of Fig. 18 we find that adding the MegaZ-LRG samples to 
the SDSS LRG and SDSS Main data considerably shrinks the 
range of biases, e.g. by more than a factor three (seven) in the 
case of erg when sampling from the lcr (2cr) confidence vol- 
ume. In combination with the SDSS Main L3 and L4 samples, 
the intervals decrease in size by roughly another 30 - 50 % (for 



the 2<r confidence volume), reaching values which are below the 
lcr statistical errors. The reduction of intrinsic alignment bias 
to subdominant levels is also evident from the comparison with 
the lcr confidence regions for constraints on the cosmological 
parameters plotted in Fig. 18. 

Hence, under the assumptions made for this prediction, and 
provided that the mean intrinsic alignment signal were accu- 
rately known and could be subtracted from the cosmic shear 
data, the uncertainty in the knowledge about the free intrinsic 
alignment parameters in (19) and (23) would be subdominant to 
the statistical errors in a CFHTLS-like survey, given the intrin- 
sic alignment constraints from the joint fit to all galaxy samples 
considered in this work. 

It is also evident from Fig. 18 that the mean bias on the pa- 
rameters Q m , erg, and wq due to our best-fit intrinsic alignment 
model is in each case appreciably smaller than the lcr statisti- 
cal errors. Consequently, it is possible that cosmology would not 
be significantly biased even if intrinsic alignments were simply 
ignored. One may ask if any of the existing cosmic shear anal- 
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Table 6. Range of possible cosmological parameter biases due to intrinsic alignments for a CFHTLS-like cosmic shear survey. 



parameter 


C stat 


sample set (1) 


sample set (2) 


sample set (3) 


sample set (4) 






lo- 


2cr 


lcr 


2o- 


la 


lo- 


lcr lcr 


n m 


0.06 


0.10 


0.44 


0.10 


0.37 


0.05 


om 


0.03 0.05 




0.11 


0.11 


0.53 


0.11 


0.42 


0.03 


0.06 


0.02 0.04 


w 


0.62 


0.49 


2.84 


0.46 


2.40 


0.40 


0.75 


0.24 0.37 



Notes. We have listed the lcr statistical error cr stat , resulting from the Fisher matrix analysis after marginalising over all remaining parameters, in 
the second column, as well as the range of biases we obtained by sampling from the lcr and 2cr confidence regions in the parameter space spanned 
by (A,/?, Mother! • In the third to sixth columns results from the fits to four sets of galaxy samples, which are also shown in Fig. 18, are given. These 
sets are (1) the combined SDSS LRG samples, (2) SDSS LRG and SDSS Main L4 and L3 samples combined, (3) SDSS LRG and MegaZ-LRG 
samples combined, and (4) SDSS LRG, MegaZ-LRG, and the SDSS Main L4 and L3 samples combined. The range of biases is defined as the 
interval containing 99 % of the distribution of biases (29). 



yses would be affected more seriously if subject to an intrinsic 
alignment signal that follows our best-fit model. Hitherto, weak 
lensing surveys have not been used in combination with photo- 
metric redshifts of individual galaxies to perform cosmic shear 
tomography, with the exception of the space-based COSMOS 
survey (Massey et al. 2007b; Schrabback et al. 2010). 

Non-tomographic surveys have much lower contamination 
by intrinsic alignments than tomographic studies because for 
a wide redshift distribution the probability of having a close 
galaxy pair is smaller than for an auto-correlation of narrow 
redshift bins, thereby lowering the II signal. The probability 
of having a large line-of-sight separation of galaxies in turn is 
smaller than for cross-correlations of a low- and high-redshift 
tomographic bin, thereby rendering the GI contribution less im- 
portant. Hence, we conclude that any intrinsic alignment sig- 
nal close to our best-fit model is irrelevant for existing surveys, 
unless they are significantly shallower, which places a typical 
galaxy of the survey sample at lower redshift and higher lumi- 
nosity (see e.g. Fig. C.2), so that the intrinsic alignment contam- 
ination becomes stronger (see also Kirk et al. 2010). 

The COSMOS survey is deeper than the CFHTLS-like sur- 
vey analysed above (i8i4 < 26.7, Schrabback et al. 2010), so that 
a substantial part of the cosmic shear signal stems from high 
redshifts. In addition, red galaxies are likely to be less lumi- 
nous on average, both effects tending to decrease the amplitude 
of intrinsic alignments. Besides, due to the small survey area of 
COSMOS, cosmological constraints are modest, further dimin- 
ishing the importance of an intrinsic alignment bias. Indeed, the 
fully tomographic analysis by Schrabback et al. (2010) did not 
detect any effects due to intrinsic alignments by excluding auto- 
correlations of tomographic bins and bright red galaxies. 

7. Conclusions 

In this work we studied intrinsic alignments in the MegaZ- 
LRG galaxy sample, investigating for the first time an early- 
type galaxy sample at intermediate redshifts up to z ~ 0.6, 
and for which only photometric redshift information is avail- 
able. We presented correlations between galaxy number densi- 
ties and galaxy shapes as a function of the transverse comov- 
ing separation of the galaxy pairs for MegaZ-LRG and two sub- 
samples at high and low redshift, as well as for several spec- 
troscopic SDSS LRG and SDSS Main samples. In combination, 
these samples comprise wide ranges in redshift (z < 0.7) and lu- 
minosity (4 mag) which have not been covered in a joint analysis 
before. 

We developed the formalism to incorporate photometric red- 
shift scatter into the modelling of the correlation function, taking 



special care of the large line-of-sight spread of physical correla- 
tions and the effect of contributions by galaxy-galaxy lensing 
and lensing magnification-shear cross-correlations, which is in- 
troduced by photometric redshift uncertainty. Our model repro- 
duces to good accuracy the scaling of the MegaZ-LRG data with 
the maximum line-of-sight separation included in the computa- 
tion of the correlation functions, as well as the relative contribu- 
tion by galaxy-galaxy lensing when varying this maximum sep- 
aration. This supports the validity of our modelling ansatz and 
justifies the use of the photometric versus spectroscopic relation 
obtained from the 2SLAQ survey to quantify the photometric 
redshift scatter. 

Moreover we discussed a correction to the redshift depen- 
dence of the widely used linear alignment model (Hirata & 
Seljak 2004, 2010). We then fitted the nonlinear version of this 
corrected linear alignment (NLA) model to all samples with a 
free overall amplitude. To allow for the transition from galaxy- 
intrinsic to matter-intrinsic correlations, we also determined a 
linear galaxy bias from galaxy clustering signals. Due to the as- 
sumption of linear biasing and the expected breakdown of the 
NLA model on small scales, we limited the analysis to comov- 
ing transverse separations r p > 6 /r'Mpc. We found that all sam- 
ples are consistent with the scaling with r p that is inherent to the 
NLA model (which is identical to the scale dependence of the 
matter power spectrum), suggesting that the alignment of early- 
type galaxies is indeed determined by the local tidal gravitational 
field. We did not test other theories of intrinsic alignments since 
the linear alignment model is, at least on the largest scales, phys- 
ically motivated, widespread, and reasonable for elliptical galax- 
ies whose shapes are expected to be subjected to tidal stretching 
by the surrounding gravitational field. 

The amplitudes of the intrinsic alignment signals that were 
obtained from the fits vary widely by more than an order of 
magnitude and are thus inconsistent with a one-parameter NLA 
model. We introduced additional power-law terms into the NLA 
model to account for an extra redshift and luminosity depen- 
dence, using in each case the power-law index as a further free 
parameter. With this three-parameter model, we demonstrated 
that the MegaZ-LRG, SDSS LRG, and SDSS Main samples un- 
der consideration are fully consistent with each other, but we 
add the caveat that the error bars for each sample are still quite 
large. We would particularly benefit from constraints from a low- 
luminosity sample that probes a larger volume since the statistics 
for the SDSS Main sample are relatively poor and thus the bright 
galaxies from the LRG samples dominate the fits. The joint fit 
to all ten samples strongly suggests an approximately linear in- 
crease of the intrinsic alignment amplitude with galaxy luminos- 
ity and is consistent with no extra redshift evolution beyond the 
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(corrected) NLA model. Adding in the new MegaZ-LRG data is 
particularly beneficial in tightening the constraints on the red- 
shift evolution of the intrinsic alignment signal due to the higher 
redshift of the samples. The normalisation of the NLA model is 
given by 0.077 p^ 1 (combining the parameters C\ and A), with a 
lcr uncertainty of roughly 10 %. 

In the joint analysis of galaxy samples, special attention 
was payed to homogenising the samples as regards the deter- 
mination of rest-frame magnitudes/luminosities and the range of 
galaxy colours. As a consequence, we re-defined the colour sep- 
arator of the red SDSS Main LA and L3 samples compared to 
Hirata et al. (2007) to avoid a leakage by blue-cloud galaxies. 
Furthermore we discarded for the joint fits about 30 % of the 
MegaZ-LRG galaxies with colours 2 (g - i) < 1.7, the resulting 
redder subsample producing slightly higher intrinsic alignment 
amplitudes. We also discussed the dependence on the employed 
k + e-corrections, affecting our results via the computation of 
luminosities and rest-frame colour cuts. 

A range of systematics tests was applied to all new and re- 
defined galaxy samples, finding no evidence for systematic ef- 
fects in the correlation functions. Residual sources of uncer- 
tainty, e.g. in the calibration of galaxy shears or due to the sta- 
tistical error on the galaxy bias, should be clearly subdominant 
to the uncertainty originating from the statistical errors on the 
correlation function measurements that were propagated into the 
errors on the intrinsic alignment parameters. 

In the linear alignment picture, the normalisation of the in- 
trinsic alignment signal is determined by the response of the in- 
trinsic shape of a galaxy to the gravitational potential of that 
galaxy at the time of its formation. Interpreting our constraints 
on intrinsic alignment parameters in this framework, we ob- 
tained no evidence for an extra redshift evolution and hence a 
time dependence of the coupling between intrinsic galaxy shape 
and gravitational potential. The dependence on luminosity we 
found can be interpreted as an increase of this coupling with 
galaxy mass. Comparing our results with Mandelbaum et al. 
(2010) who analysed blue galaxy samples out to similar red- 
shifts, it is evident that intrinsic alignments depend strongly on 
the galaxy type. Whether there is a clear dichotomy between 
early- and late-type galaxies or a more continuous transition with 
galaxy colour, which the comparison between the MegaZ-LRG 
samples with and without colour cut may hint at, is still to be 
investigated. 

In a Fisher matrix analysis, we predicted the bias on cosmo- 
logical parameters that results from our best-fit intrinsic align- 
ment model. To this end, we used sets of luminosity functions 
measured by Faber et al. (2007) in order to derive the frac- 
tion of early-type galaxies and their luminosity distribution at 
each redshift for a given apparent magnitude limit of the survey. 
Assuming zero intrinsic alignments for blue galaxies (without 
uncertainty), we then computed the expected intrinsic-ellipticity 
(II) and shear-intrinsic (GI) correlations, sampling the parame- 
ters of the NLA model from the confidence volume of our in- 
trinsic alignment fits. The accuracy of this approach is limited 
by the substantial extrapolation of the luminosity function data 
to faint and high-redshift galaxies, as well as the strong scatter in 
luminosity function parameters obtained from different works. 

For a fully tomographic CFHTLS-like survey both the mean 
bias and the scatter in the bias due to the uncertainty in the intrin- 
sic alignment model are smaller than the predicted lcr cosmo- 
logical parameter errors, and we similarly expect subdominant 
systematic effects by intrinsic alignments for other cosmic shear 
surveys performed hitherto. The MegaZ-LRG data were crucial 
in reducing this systematic uncertainty. However, if external pri- 



ors on cosmological parameters, e.g. from the cosmic microwave 
background, are employed, which is likely to be the case in prac- 
tice, the significance of the bias due to intrinsic alignments will 
be higher. 

For future ambitious weak lensing surveys such as Euclid, 
which has roughly comparable depth to CFHTLS but consid- 
erably higher statistical power (Laureijs et al. 2009), the same 
intrinsic alignment signal would constitute a severe systematic, 
and marginalising over the uncertainty in intrinsic alignment 
parameters would significantly degrade constraints on cosmol- 
ogy. However, the gradually increasing precision requirements 
by planned cosmic shear surveys are likely to be matched by in- 
trinsic alignment studies that continuously improve and consoli- 
date the models. Hence, constraints on intrinsic alignment mod- 
els as provided by this work and succeeding ones will prove most 
useful, for instance to define prior ranges in nuisance parametri- 
sations of intrinsic alignment signals (e.g. Bridle & King 2007). 

A straightforward test of the intrinsic alignment model ob- 
tained in this work could be provided by a similar analysis of 
galaxy shape correlations from the same set of samples. This 
way one can verify whether the II signal agrees with the pre- 
diction by the linear alignment paradigm (based on the GI sig- 
nal in this work), and whether the extra redshift and luminosity 
dependencies of the II signal are consistent with the present re- 
sults. Furthermore, to obtain unbiased intrinsic alignment mea- 
surements for a wide range of galaxy colours is another impor- 
tant but challenging task because many central selection criteria 
such as shape measurement quality, photometric redshift scat- 
ter, or spectroscopic redshift failure rates depend strongly on the 
galaxy type. 

Measurements of the type presented in this paper are re- 
stricted to quasi-linear scales although both cosmic shear and 
galaxy evolution studies have a vital interest in intrinsic align- 
ments in the deeply nonlinear regime. A possible way forward 
would be the usage of a halo model approach for both galaxy 
bias and intrinsic alignment signals (see e.g. Schneider & Bridle 
2010). However, note that, similar to the galaxy bias, observa- 
tional constraints on intrinsic alignments may ultimately be lim- 
ited by an intrinsic scatter that cannot be accounted for by a de- 
terministic model. 

The generalisation of intrinsic alignment measurements to 
photometric redshift data has opened up a new regime of data 
sets which could be exploited to constrain intrinsic alignment 
models. For instance, all present or upcoming cosmic shear sur- 
veys with redshift information, or at least subsamples of them 
with low photometric redshift scatter, could be suited, thereby 
automatically extending the baselines in redshift and luminos- 
ity to scales most relevant for weak lensing. The higher the 
photometric redshift scatter, the more important become galaxy- 
galaxy lensing and magnification contributions to the observed 
correlation functions, so that e.g. at some point cosmology will 
have to be varied in the intrinsic alignment analysis as well. Then 
it might be more fruitful to instead consider simultaneous mea- 
surements of galaxy shape and number density correlations in 
the manner proposed by Bernstein (2009) and investigated in 
Joachimi & Bridle (2010). 
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Appendix A: Three-dimensional number 
density-intrinsic shear correlation function 

In this appendix, we derive the three-dimensional number 
density-intrinsic shear (gl) correlation function, detail the inclu- 
sion of photometric redshift scatter into the formalism, and es- 
tablish an approximate relation between the three-dimensional 
correlation function and the angular power spectrum. 

A.1. Correlation function for exact redshifts 

We define the three-dimensional correlation function between 
the galaxy density contrast 6 g and the radial intrinsic shear yi >+ 
as 

t g i(r p ,Il,z) = (s g {o,x(z)- |,z) yi, + (*L,;Kz)+ f ,z)) ,(A.i) 

for a given mean redshift z of the galaxy pairs correlated. Here 
we introduced a three-dimensional comoving separation vector 
x which has a line-of-sight component II = x\\. Its transverse 
components are denoted by x±, with modulus r p = \x±\. The 
first argument of both 6 g and yi j+ denotes the position on the 
sky, the second the position along the line of sight, and the third 
quantifies the epoch, given in terms of the redshift. Note that a 



line-of-sight separation LI + implies that 6 g and are not 
measured at precisely the same epoch, contrary to what we have 
written in (A.l). However, as n is small compared to the co- 
moving distance x(z) to me galaxies under consideration, this 
approximation holds to good accuracy. 

Following Hirata & Seljak (2004), the radial component of 
the intrinsic shear is measured with respect to x±, and without 
loss of generality we can choose the coordinate system such that 
Ti,+ - 7i,i • Note that in the majority of weak lensing studies y + is 
defined as the tangential component of the shear. Measuring ra- 
dial instead of tangential alignment implies a change of sign, so 
that e.g. the galaxy-galaxy lensing signal which we consider is 
negative. Denoting Fourier variables by a tilde, one can construct 
in analogy to the matter density contrast a three-dimensional in- 
trinsic convergence <5i 12 via 

Hk) = e 21 ^ Uk) , (A.2) 

where ip is the polar angle of k ± , i.e. the projection of the wave 
vector onto the plane of the sky. We will denote the line-of-sight 
component of k by k\\. 

Then one can write the correlation function by Fourier trans- 
forming (A.l) as 

&<r„IM = J f ^ I e-'« '(};(t'. I )f w (t. z ))(A.3) 

Inserting the definition of the three-dimensional gl power spec- 
trum, 

(6;(k',z) fi(*,z)) = (2tt) 3 dg) (k - k') P gI (k,z) , (A.4) 
and subsequently integrating (A. 3) over k' yields 

n, z) = f e ~ lkx ^ *) < A - 5 > 

where in order to arrive at the third equality, the definition of 
the second-order Bessel function of the first kind was used. In 
this derivation it was implicitly assumed that the intrinsic shear 
field does not feature B-modes, as is for instance the case for the 
linear alignment paradigm. 

One can now integrate over the line of sight, making use of 
the definition of the Dirac delta-distribution, to obtain the pro- 
jected gl correlation function as employed by Mandelbaum et al. 
(2006), Hirata et al. (2007), and Mandelbaum et al. (2010), 

dnf gl (r p ,Yl,z) (A.6) 

CO 

r°° dk ± k ± 

= - —z Ji(kj_r p ) P s i(k ± ,z) ■ 

J0 2 7T 

12 The intrinsic convergence Si can best be interpreted as the conver- 
gence of the intrinsic shear field for infinitesimal slices in comoving 
distance. Intrinsic shear field, projected intrinsic convergence, and Si 
have been defined in exact analogy to the cosmic shear formalism. For 
details we refer the reader to Joachimi & Bridle (2010) who employ the 
same notation except for using K\ instead of Si. 



30 



B. Joachimi et al,: Constraints on intrinsic alignment contamination of weak lensing surveys using the MegaZ-LRG sample 



Real data cannot provide the correlation function for arbitrar- 
ily large line-of-sight separations, so that a truncation of the in- 
tegral in (A. 6) is necessary. This formula is still applicable if 
one can stack observations for all values of n for which galaxy 
pairs carry a signal. While this can easily be achieved for spec- 
troscopic observations, photometric redshift scatter smears the 
signal in II such that a cut-off n max needs to be taken into ac- 
count explicitly in the modelling. Of course it would be possible 
to compute the observed correlations out to very large n max , but 
this way many uncorrelated galaxy pairs would enter the corre- 
lation function, thereby decreasing the signal-to-noise dramati- 
cally. 

Instead, we proceed from (A. 5) by assuming that £ g i is a real 
function, and write 



£ g i(r p ,n )Z ) = 
Jo n Jo 



(A.7) 



dk ± k ± 
2n 



Jiik^p) pJyjkT+k\,z \ cos(fc||II) . 



As can be seen from this equation, £ g i is an even function in both 
r p and n, so that it is sufficient to compute just one quadrant. 
Note that by definition r p > 0, whereas II can also attain negative 
values. 

Equation (A.7) yields the three-dimensional gl correlation 
function for exact or, to good approximation, spectroscopic red- 
shifts. For the model described in Sect. 4 with SuperCOSMOS 
normalisation and b g - 1, we plot ^ g i(r p ,H, z) for z ~ 0.5 in 
Fig. A.l, bottom panel. As expected, the correlation is strongest 
for small separation, in particular for |ll| close to zero. If spec- 
troscopic data is available, essentially all information is captured 
when a cut-off n max = 60/z'Mpc is used in the integration 
(A. 6), as e.g. in Mandelbaum et al. (2010). Due to the defini- 
tion (A.l), the gl correlation function measures the radial align- 
ment of the galaxy shape with respect to the separation vector 
of the galaxy pair considered. Therefore the correlation func- 
tion vanishes for all II at r p = since then the separation vec- 
tor points along the of sight. Note that the contours do not ap- 
proach the II = 0-axis asymptotically, but cross this line at some 
value of r p , as expected for a differentiable correlation function. 
Throughout these considerations we have not taken the effect of 
redshift-space distortions into account. 
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Fig. A.l. Three-dimensional gl correlation function as a function of 
comoving line-of-sight separation If and comoving transverse sepa- 



ration r„ at z 



0.5. 



Contours are logarithmically spaced between 
10~ 2 (yellow) and 10~ 6 (black) with three lines per decade. Top panel: 
Applying a Gaussian photometric redshift scatter of width 0.02. Bottom 
panel: Assuming exact redshifts. Note the largely different scaling of 
the ordinate axes. The galaxy bias has been set to unity, and (6) with 
SuperCOSMOS normalisation has been used to model Pgi in both cases. 
Redshift-space distortions have not been taken into account. 



A.2. Incorporating photometric redshifts 

Photometric redshift errors cause the observed correlation func- 
tion to be a 'smeared' version of (A.7), introducing a spread es- 
pecially along the line of sight but to a lesser extent also in trans- 
verse separation (because an uncertain redshift is used to convert 
angular separation to physical separation). If we denote quanti- 
ties determined via photometric redshifts by a bar, the actually 
measured three-dimensional correlation function reads 



(A.8) 



J"dz m j dr p J dll p(r p ,n,z m \r p ,Il,z m ) ^ g i(r p ,Tl,z m ) , 

where z m denotes the mean redshift of the galaxy samples used 
for the number density and the shape measurement. Here, p 
is the probability distribution of the true values of r p , II, and 
z m , given photometric redshift estimates of these quantities. In 
words, (A.8) means that in order to obtain the observed correla- 
tion function, we integrate over £ g i as given in (A.7), weighted 
by the probability that the true values for separations and red- 
shift actually correspond to the estimates based on photometric 
redshifts. 



The direct observables for this measurement are the red- 
shifts of the two galaxy samples under consideration, z\ and zi, 
and their angular separation 0. The sets of variables {z\,zi,0) 
and (r p , II, z m ) are related via a bijective transformation. Writing 
(A.8) in terms of the other set of variables, one obtains 



C'^n,z m ) 



(A.9) 



j dzi J dz 2 J d0p(z 1 ,Z2,O\z 1 {z m ,Tl},z2 [z m ,Tl} ,6[z m ,r p }) 

X £gl (r p {Zl , Z2, 0} , 11 {Zi , Z 2 ) , Zm {Zl , z 2 }) 
= J dzi J dz 2 Pn {z\\l\ [z m , A}) p t (z2\Z2 {im, A}) 

X £gl { r P {^1 ' 22, 0(Zm, rp)} , II {Z\ ,Z2},Z m {Zl , Z2}) ■ 

In the second step it was assumed that the probability distribu- 
tions of zi, Z2, and are mutually independent, and that 6 is ex- 
actly known, i.e. p(0\0) = Sd(6 - 0)- We have introduced differ- 
ent redshift probability distributions for the galaxy sample with 
number density information p„ and the one with shape informa- 
tion p e . All quantities related to photometric redshifts have been 
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expressed in terms of the arguments of the correlation function 
on the left-hand-side. 

We make use of the following approximate relations between 
the two triples of variables, 



be truncated for reasons of a good signal-to-noise ratio, so that 
in the case of photometric redshifts part of the signal is lost. 
Therefore it is crucial to repeat the same steps applied to the 
data also to the model and use the same cut-off U max in (A. 6). 



Zm = 2 fa + zi) ; 

r n ~ Ox(Zm) ; 
c 



•p 
n 



H(z m ) 



(Z2-Z\) , 



where H(z) is the Hubble parameter. Note that the same trans- 
formations have been used to bin the observational data in terms 
of redshift, transverse and line-of-sight separation. With this 
equation for II, in combination with the assignment of prob- 
ability distributions in (A.9), we have introduced the conven- 
tion that II > means that the galaxy from the density sample 
is at lower redshift than the galaxy from the shape sample. If 
and only if the distributions for the density and the shape sam- 
ple are identical, which we assume throughout this work, the 
correlation function remains symmetric with respect to II, i.e. 

^,n,z ffl )^; o, (r P ,-n ; z m ). 

With these equations at hand, one can also write down the 
inverse transformation of (A. 10), which is needed to evaluate 
(A.9), 



(A.10) ^ 3 R e i a tion to angular power spectra 

We now derive a relation between the three-dimensional gl cor- 
relation function in the presence of photometric redshift scatter 
and the angular power spectrum, which proves most convenient 
to compute in practice. Inserting (A.5) into (A. 12), one can 
write 



^ ot (^n,z m ) = 



- I dzi I dz 2 Pn I Zl I Zm — | Pe 



2c 



Zl I Zm + 
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where we employed z m = (z\ + Zi)l2 as a shorthand notation. 
Making use of 6 = r p /x(Zm), see (A. 11), and defining the angular 
frequency £ = k ± x (z m ). one obtains 
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Then (A.9) can be expressed as 
^ 0t (r„n,Zm) = 
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We then transform the integration variables {zi,Z2} to 
{z m , Az = Z2 - Zi). Note that the determinant of the Jacobian of 
this transformation is unity. We apply Limber's approximation, 
which in this case can be written as 



Note that the absolute value for Z2 - Zi has been introduced in 
the second argument of £ g i, which is possible since it is an even 
function in this argument. The integrals in (A. 12) run over the 
full range of spectroscopic (exact) redshifts. As a consequence, 
|Z2 - Zi| in the second argument of £ g i can obtain relatively large 
values, leading to very large II » lOO/r'Mpc. However, the 
spectroscopic £ g i becomes very small for large n, so that the 
integrand in (A. 12) can safely be set to zero in this case. 

Still, any sizeable photometric redshift scatter leads to a con- 
siderable spread of the three-dimensional correlation function in 
n, as can be seen in Fig. A.l. Assuming a Gaussian photomet- 
ric redshift scatter with width 0.02 around every true redshift, 
the strong signal concentrated at small n and r p < lO/z^Mpc 
in the spectroscopic case is scattered along the line of sight, so 
that the values of ^ ot at n > 200/T 1 Mpc are still more than a 
per cent of those at n = for any r p . In contrast, we find that 
the net scatter of signal between different transverse separations 
is negligible. Hence, in principle the projected correlation func- 
tion (A. 6) does not change when using photometric instead of 
spectroscopic redshift information as long as the complete range 
of n for which a signal is measured enters the line-of-sight in- 
tegration. However, in practice the line-of-sight integral has to 
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Here we have assumed that the two redshift probability distri- 
butions are sufficiently broad and have similar forms, so that an 
evaluation at z m instead of z m ± Az/2 does not change the results 
significantly. Since the photometric redshifts on which the dis- 
tributions are conditional encapsulate the dependence of ^j 10 ' on 
the line-of-sight separation n, we do not extend this approxima- 
tion to the second argument. Equation (A.14) thereby simplifies 
to 
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where in order to arrive at the second equality, we integrated over 
Az. The resulting Dirac delta-distribution renders the k\\ integra- 
tion trivial. Making use of the expressions dz m = dxH (z m ) /c 
and p(z) = p(x) d^/dz, one obtains the result 
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AT 



where we wrote the normalisation in the notation of Hirata & 
Seljak (2004) in which C\ is an arbitrary constant. The par- 
tial derivatives are with respect to comoving coordinates, and 
<S> p (x) = cD(;c,Zp) denotes the 'primordial' potential, i.e. the lin- 
ear gravitational potential evaluated at the epoch of galaxy for- 
mation, at a redshift z p well within the matter-dominated era. For 
ease of notation we have omitted a smoothing of the gravitational 
potential on galactic scales in (B.l) which can be implemented 
by a simple cut-off of high wavenumbers in Fourier space, see 
Hirata & Seljak (2004). These authors used the relations (B.l) in 
their derivation of the intrinsic alignment power spectra for the 
linear alignment model, which we closely follow. 

In a first step, the primordial gravitational potential is related 
to the matter density contrast via the Poisson equation 



(B.2) 



where in the last step we implicitly defined the projected gi 
auto-correlation power spectrum C g i. In addition to the angu- 
lar frequency, we have written the photometric redshifts z\ and 
Z2, which characterise the redshift distributions entering C g i, ex- 
plicitly as arguments. Note that Limber equations, such as (5), 
in general hold only approximately, the range of validity being 
the more limited the narrower the kernels in the line-of-sight in- 
tegration (e.g. Simon 2007). 

We have verified that the calculations of the three- 
dimensional gi correlation function according to (A. 12) and 
(A. 17) agree within the numerical accuracy. The latter can 
be computed much more efficiently by computing the angular 
power spectrum via Limber's equation and then using Hankel 
transformations to obtain the correlation function ^ ot (r p , ft, z m ) 
via (A. 17), employing the transformation (A. 10). One can pro- 
ceed likewise to obtain analogous expressions for the gg signal. 
Galaxy-galaxy lensing vanishes if the density field probed by the 
galaxy distribution and the source galaxies on whose images the 
gravitational shear is measured are located at exactly the same 
redshift. Thus one cannot proceed with the same formalism as 
used to derive the gi contribution (see the assumptions under- 
lying the definition (A.l)), but must instead incorporate redshift 
probability distributions from the start, again arriving at an ex- 
pression analogous to (A. 17). 

Appendix B: Redshift dependence of the linear 
alignment model 

In this appendix we re-derive the redshift dependence of the lin- 
ear alignment model, obtaining a different result than Hirata & 
Seljak (2004), but being in full agreement with Hirata & Seljak 
(2010). Practically all attempts at constructing a physical de- 
scription for intrinsic alignments are based on the linear align- 
ment model originally suggested by Catelan et al. (2001). They 
assumed that the shape of the luminous part of a galaxy exactly 
follows the shape of its host halo, and that the ellipticity of the 
latter is determined by the local tidal gravitational field of the 
large-scale structure. 

The simplest possible form allowed by the assumptions 
made above is a linear relation between the intrinsic shear and 
the gravitational field, given by (Catelan et al. 2001) 
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= 47iG Pm (z)(l+z)-' 6(x,z P ) , 



where V 2 denotes the comoving Laplacian. This expression is 
Fourier-transformed, yielding 



<&(*,z) = -47rGp m (z) (1 +z)~ 2 r 2 6(k,z) 



(B.3) 



The growth factor D(z) quantifies the dependence of the mat- 
ter density contrast on redshift in the limit of linear structure 
formation, and is normalised to D(z) = (1 + z) _1 during matter 
domination (Hirata & Seljak 2004) 13 . Restricting (B.3) to linear 
scales, one obtains the ratio 



d>(*,z) (l+z)D(z) 



<D(ft,z p ) (l+Zp)£>(z P ) 



(l+z)D(z), 



(B.4) 



where in the last step we made use of the fact that z p lies in the 
matter-dominated era. Inserting (B.4) into (B.3), and considering 
the linear regime, one arrives at the relation between primordial 
potential and linear matter density contrast, 



%(k) = - 



4nG p m (z) 
D(z) (1 + z) 3 



k S hn (k,z) . 



(B.5) 



This expression differs from the result given in Hirata & Seljak 
(2004), eq. (14), by an additional factor (1 + z)~ 2 . This discrep- 
ancy was also found in other re-derivations (Hirata & Seljak 
2010; R. Bean, I. Laszlo; priv. comm.). 

Hirata & Seljak (2004) inserted their eq. (14) into (B.l) and 
then computed the three-dimensional intrinsic shear (II) power 
spectrum and the intrinsic shear-matter cross-power spectrum. 
Neglecting source clustering but otherwise following their steps 
in exact analogy, we obtain 
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On 



Psi(k,z) = -C\ p CT — — PsM(k,z) , 
D(z) 



(B.6) 
(B.7) 



where in this work we use the full nonlinear matter power spec- 
trum on the right-hand side instead of the linear power spec- 
trum, as written here in the original form of the linear alignment 
model. If Hirata & Seljak (2004), eq. (14), were employed in- 
stead of (B.5), (B.6) would have an additional term (1 + z) 4 , and 



13 Note that this normalisation of the growth factor differs from the 
one used in the remainder of this work. However, this difference only 
affects the amplitude C\ which is an arbitrary constant. Our default nor- 
malisation of the growth factor is in line with Bridle & King (2007) 
whose value for C\ we have adopted. 



33 



B. Joachimi et al.: Constraints on intrinsic alignment contamination of weak lensing surveys using the MegaZ-LRG sample 




Fig. C.l. Red galaxy fraction f r (top panel) and mean luminosity of red 
galaxies (bottom panel) as a function of redshift for a magnitude limit 
r lim = 25. We compare the results for the sets of luminosity functions 
provided by Faber et al. (2007, black solid curves), Brown et al. (2007, 
dark grey solid curves), Willmer et al. (2006, light grey solid curves), 
Giallongo et al. (2005, black dotted curves), and Wolf et al. (2003, grey 
dotted curves). The results for Faber et al. (2007) luminosity functions 
used with the value of B - r expected at z = 1 are indicated by the cross 
in each panel. 

(B.7) an additional term (1 +z) 2 , in the numerator. These modifi- 
cations would correspond to a shift by -2 in pother in our models 
(19) and (23). 

Appendix C: Volume density and luminosities of 
red galaxies 

To make realistic predictions for the intrinsic alignment contami- 
nation of cosmic shear surveys, we must specify, at each redshift, 
the distribution of galaxy luminosities that enter (19). Since this 
intrinsic alignment model only holds for red galaxies, we addi- 
tionally must estimate the fraction of early-type galaxies in the 
total weak lensing population as a function of redshift. In this pa- 
per, both quantities are determined using fits to the observed lu- 
minosity functions given in Faber et al. (2007). In this appendix, 
we present technical details about these calculations, assess the 
sensitivity of our results to this particular luminosity function 
choice, and provide data that can be used to forecast the intrin- 
sic alignment contamination of other cosmic shear surveys (with 
different limiting magnitudes) besides that discussed in the main 
text. In all cases, we are extrapolating the luminosity functions 
to fainter magnitudes at a given redshift relative to the samples 
used to determine the luminosity function. 

We employ the Schechter luminosity function parameters for 
red galaxies from Faber et al. (2007), where <p* and M* are given 
as a function of redshift, and where the faint-end slope is fixed at 



a = -0.5. While we consistently use magnitudes in the r band, 
Faber et al. (2007) determine M* in the B band. Therefore we 
the estimate rest-frame B - r colour from the tables provided in 
Fukugita et al. (1995), finding B - r — 1.32 for ellipticals. This 
conversion from B to r takes into account that Faber et al. (2007) 
give B band magnitudes in the Vega-based system, whereas this 
work uses AB magnitudes throughout. Furthermore, we have 
assumed r « r', where r' is the filter listed by Fukugita et al. 
(1995). This assumption should hold to good accuracy 14 for typ- 
ical colours of the galaxies in our samples, i.e. 0.2 < r — i < 0.6. 

For early-type galaxies, B - r shows little evolution between 
z = and z ~ 1 (Bruzual & Chariot 2003), so we assume the 
rest-frame colour to be constant in this redshift range, which we 
check via the following procedure. Since the Sloan g filter cov- 
ers a similar wavelength range to the B band (although the peaks 
of the transmission curves differ, see Fukugita et al. 1995 for de- 
tails), we use the evolution of g - r as determined from the Wake 
et al. (2006) templates as an approximation for the redshift de- 
pendence of B - r. We find a shift of 0.15 mag from z = to 
Z = 1, which has significantly less effect on our results than em- 
ploying different observational results for luminosity functions, 
see Fig. C. 1 and the corresponding discussion below. Finally, we 
correct for the fact that Faber et al. (2007) have computed abso- 
lute magnitudes assuming a Hubble parameter h = 0.7 while we 
give absolute magnitudes in terms of h = 1 . 

With all of these caveats, the limiting absolute B band mag- 
nitude at redshift z for a given apparent magnitude limit in the r 
band is given by 



^min(£? Him) — Him 



1 Mpc 



+ {B-r), 
(CI) 



where k r?le d(z) is the ^-correction of red galaxies for the r band 
(Wake et al. 2006). In line with our convention for absolute mag- 
nitudes, the luminosity distance Di is computed with h = 1 . If 
absolute magnitudes are given for other values of the Hubble 
parameter, like e.g. in Faber et al. (2007), we convert these ac- 
cordingly. The limiting absolute magnitude from (C.l) can then 
be transformed into the minimum luminosity entering (27) and 
(28), 
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(C2) 



^minfc Him) 

Lo(z) 

where Mq(z) denotes the rest-frame absolute magnitude -22, 
evolution-corrected to redshift z using the redshift dependence 
of M* from Faber et al. (2007), which is given by - 1 .2z. Note 
that this dependence accounts for the redshift evolution in the B 
band, but since B - r is nearly constant as a function z, we can 
also apply the correction (to good approximation) to r band mag- 
nitudes. Denoting the luminosity corresponding to M* by L*, we 
obtain for (28) the expression 

L" (z) Uo(z)j r(a+ 1, 

where the incomplete Gamma function T(a, x) — J dyy a ~ l e~ 
was introduced. Analogously, we arrive at 



(C3) 



"V,red(z, Him) = <f>*(z) T Or + 1 , 



^min (Z, Him) 
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(C4) 



14 See http : //www . sdss . org/dr6/algorithms/ 
jeg_photometric_eq_drl .html#usno2SDSS for the transformation 
equation between r and r' . 
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Fig. C.2. Comoving volume number density of red galaxies n\ T and mean luminosity (L) /Lq of red galaxies as a function of limiting magnitude 
rii m and redshift z. Left panel: Comoving volume number density of red galaxies n\ f in units of 1CT 4 Mpc~ 3 . Contour values range between 10" 3 
in the upper left and 50 at the bottom. In the upper left corner of the panel fiy.r ~ 0. Right panel: Mean r band luminosity. Contour values range 
between 0. 16 in the lower part to 5 in the upper left corner. In both panels the red dashed line marks the values for nim = 25, the limiting magnitude 
we employ in our calculations. Decades in contour values are indicated by the black solid lines. 



for the comoving volume density of red galaxies entering (27). 

In addition to the luminosity functions from Faber et al. 
(2007), we also consider fitted Schechter parameters presented 
in Giallongo et al. (2005), as well as the sets of luminosity func- 
tions published by Wolf et al. (2003), Willmer et al. (2006), and 
Brown et al. (2007). We determine fit functions to the redshift 
dependence of both M* and (f>* for the latter three works because 
we have to extrapolate beyond the range of redshift analysed 
therein. We use linear functions for M* and various functional 
forms with two to three fit parameters for <p*, but note that since 
the fits rely on only five to six data points, the extrapolation has 
considerable uncertainty. All five references give B band lumi- 
nosity functions, but the magnitude system and the convention 
for h vary, as well as the redshift ranges covered and the defini- 
tion of red galaxies. 

In Fig. C.l the red galaxy fraction f r and the mean luminos- 
ity (L) /Lq for ri; m = 25 are plotted as a function of redshift, 
making use of the different luminosity functions. We find fair 
agreement between the results based on Faber et al. (2007) and 
Brown et al. (2007), while the mean luminosities derived from 
Willmer et al. (2006) already deviate considerably at high z al- 
though Faber et al. (2007) and Willmer et al. (2006) partly use 
the same data. The Wolf et al. (2003) luminosity functions pro- 
duce significantly lower f r and higher (L) at low redshifts which 
is caused by the very different value for the faint end slope, 
a = +0.52. We note that one of the three fields chosen by Wolf 
et al. (2003) contained two massive galaxy clusters, so that the 
large-scale structure in this field could strongly influence the lu- 
minosity function in particular of early-type galaxies. However, 
small red galaxy fractions can be compensated by higher lumi- 
nosities in (19), so that even the Wolf et al. (2003) luminosity 
functions may yield intrinsic alignment signals of similar mag- 
nitude to the results of e.g. Faber et al. (2007). 

Applying the formalism to luminosity functions from 
Giallongo et al. (2005), we obtain very high f r at low redshift, 
which is clearly inconsistent with the other observations. The 



Table C.l. Mean luminosity of red galaxies and comoving volume 
number density of red galaxies as a function of redshift and limiting 
r band magnitude Him- 
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0.4 
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0.8 
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20 


0.16 


0.29 
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0.15 
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0.26 


0.37 


0.56 


1.09 


25 


0.15 
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45.79 
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25.38 


17.66 


11.57 


5.84 


27 


58.02 
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20.36 


14.63 


9.21 


28 
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46.81 
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28.80 
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16.70 


11.73 



Notes. In the upper section values for (L) /Lq are given; in the lower sec- 
tion «v.r is given in units of 10~ 4 Mpc~ 3 . This data is plotted in Fig. C.2. 



red galaxy sample used for the fits of Giallongo et al. (2005) is 
very small and contains only galaxies with z > 0.4. While the re- 
sulting fit function captures the pronounced decrease in number 
density for high redshifts that Giallongo et al. (2005) observe, 
it can obviously not be used at z < 0.4. In conclusion, we find 
that the sets of luminosity functions by Faber et al. (2007) who 
jointly analyse galaxy samples from four different surveys pro- 
duce reasonable red galaxy fractions and luminosities although 
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the uncertainty in the values of f r and (L) at any given redshift 
is still large. 

In Fig. C.2, the comoving volume density of red galaxies n v ,r 
and the mean luminosity of red galaxies in terms of Lq are plot- 
ted as a function of both ri; m and redshift, using the default set 
of luminosity functions from Faber et al. (2007). For a fixed 
magnitude limit, (L) increases with redshift while « v ,r decreases 
strongly at high redshifts. Both tendencies are more pronounced 
if the apparent magnitude limit is brighter. At low redshift, e.g. 
for z < 0.3 at ri; m = 25, a large number of faint blue galax- 
ies are above the magnitude limit and cause f r to diminish for 
z approaching zero (see Fig. C.l) although n v ,r continues to in- 
crease. This behaviour might change when explicitly taking into 
account the size cuts inherent to weak lensing surveys, but in any 
case, galaxies at these low redshifts constitute only a small frac- 
tion of the total survey volume and are expected to have a low 
luminosity and hence low intrinsic alignment signal on average, 
so they are unlikely to affect our results severely. 

The data shown in Fig. C.2, with the corresponding numbers 
collected in Table C.l, can be used in combination with the in- 
trinsic alignment model fits presented in Sect. 5.5 to predict the 
effect of intrinsic alignments on cosmic shear surveys. The mean 
luminosity as a function of redshift for a given magnitude limit 
Him can be inserted into the luminosity term in (19), which con- 
stitutes a fair approximation as long as values of /3 close to unity 
(which includes our best-fit value /3 = 1 . 1 3) are probed. Together 
with an overall redshift distribution ptotiz) for the cosmic shear 
survey under consideration, ny,r can be read off and used with 
(27) to compute the fraction of red galaxies as a function of 
redshift. Again, we emphasise the limitations of this approach 
which relies on a substantial amount of extrapolation in lumi- 
nosity, especially for fainter limiting magnitudes, and which is 
subject to the large intrinsic uncertainty in the different sets of 
luminosity functions. 
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